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Applications of perturbative QCD to deeply virtual Compton scattering and hard 
exclusive elcctroproduction processes require a generalization of the usual parton 
distributions for the case when long-distance information is accumulated in non- 
diagonal matrix elements of quark and gluon light-cone operators. I describe two 
types of nonpcrturbativc functions parametrizing such matrix elements: double 
distributions and skewed parton distributions. I discuss their general properties, 
relation to the usual parton densities and form factors, evolution equations for 
both types of generalized parton distributions (GPD), models for GPDs and their 
applications in virtual and real Compton scattering. 
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1 Introduction 



The standard feature of applications of pcrturbative QCD to hard processes 
is the introduction of phcnomcnological functions accumulating information 
about nonperturbative long-distance dynamics. The well-known examples 
are the parton distribution functional f p /H{ x ) used in perturbative QCD 
approaches to hard inclusive processes, and distribution amplitudes ^{x), 
ippf(xx,X2,X3) which naturally emerge in the asymptotic QCD analyses 2_7 of 
hard exclusive processes. More recently, it was argued that the gluon dis- 
tribution function f g (x) used for description of hard inclusive processes also 
determines the amplitudes of hard exclusive. JJip (Ref. 8) and p- meson (Ref. 9) 
electroproduction. Later, it was proposedlilrEl to use another exclusive process 
of deeply virtual Compton scattering j*p — > 77/ (DVCS) for measuring quark 
distribution functions inaccessible in inclusive measurements (earlier discus- 
sions of nonforward Compton-likc amplitudes 7*p — * 7*p' with a virtual pho- 
ton or Z° in the final state-can be found in Refs. 12-14). The important 
feature (noticed long ago E3tj) is that kinematics of hard elastic electropro- 
duction processes (DVCS can be also treated as one of them) requires the 
presence of the longitudinal (or, more precisely, light-cone "plus" ) component 
in the momentum transfer r = p — p' from the initial hadron to the final: 
r+ = Qp + . For DVCS and p-electroproduction in the region Q 2 ^> \t\,m 2 H , the 
longitudinal momentum asymmetry (or "skewedness" ) parameter £ coincidesll3 
with the Bjorken variable xgj = Q 2 /2(pq) associated with the virtual photon 
momentum q. This means that kinematics of the nonperturbative matrix ele- 
ment (p'\ ■ ■ ■ \p) is asymmetric (skewed). In particular, the gluon distribution 
which appears in hard elastic diffraction amplitudes differs from that studied 
in inclusive processesEJjO In the latter case, one has a symmetric situation 
when the same momentum p appears in both brackets of the hadron matrix 
element MjjL. ■ ■ \p). Studying the DVCS process, one deals with essentially off- 
forwardt^EH or nonforwardxJi't3 kinematics for the matrix element {p'\ . . . \p). 
Perturbative quantum chromodynamics (PQCD) provides an appropriate the- 
oretical framwork. The basics of the PQCD approaches incorporating the new 
generalized parton distributions (GPDs) were formulated in Refs. 10,11,16-19. 
A detailed analysis of PQCD factorization for hard meson electroproduction 
processes was given in Ref. 20. 

Our goal in the present paper is to review the formalism of the generalized 
parton distributions based on the approach outlined in our papers 16-18,21-25. 

Its main idea is that constructing a consistent PQCD approach for ampli- 
tudes of hard exclusive electroproduction processes one should treat the initial 
momentum p and the momentum transfer r on equal footing by introducing 
double distributions (DDs) F(x,y), which specify the fractions of p + and r + 
carried by the active parton of the parent hadron. These distributions have 
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hybrid properties: they look like distribution functions with respect to x and 
like distribution amplitudes with respect to y. 

The other possibility is to treat the proportionality coefficient ( = r + /p + 
as an independent parameter and i ntro duce an alternative description in terms 
of nonforward parton distributionst-?Mi (NFPDs) J-q{X; t) with X = x + y( be- 
ing the total fraction of the initial hadron momentum p + taken by the initial 
parton. The shape of NFPDs explicitly depends on the parameter £ char- 
acterizing the skewedness of the relevant nonforward— matrix element. This 
parametrization is similar to that proposed by X. Jil^3'tlHl3 who introduced 
off-forward parton distributions (OFPDs) H(x,£,;t) in which the parton mo- 
menta and the skewedness parameter £ = r + /2P + are measured in units of the 
average hadron momentum P = (p+p')/2. OFPDs and NFPDs can be treated 
as particular forms of skewed parton distributions (SPDs). One can also in- 
troduce the version of DDs ("a-DDs", see Ref. 22) in which the active parton 
momentum is written in terms of symmetric variables: k + = xP + + (l+a)r + /2. 

The paper is organized as follows. In Sec. 2, I recall the basic properties 
of "old" parton distributions, i.e., I discuss the usual parton densities f(x) and 
the meson distributions amplitudes <p{a). In Sec. 3, I consider deeply virtual 
Compton scattering as a characteristic process involving nonforward matrix 
elements of light-cone operators. I introduce double distributions f(x, a; t) and 
discuss their general properties. The alternative description in terms of skewed 
parton distributions H(x,£;t) is described in Sec. 4. Models for double and 
skewed distributions based on relations between GPDs and the usual parton 
densities are constructed in Sec. 5. The evolution of GPDs at the leading 
logarithm level is studied in Sec. 6. In Sec. 7, I discuss recent studies of 
DVCS amplitude at twist-2 and twist-3 level. In Sec. 8, the GPD formalism 
is applied to real Compton scattering at large momentum transfer. In the 
concluding section (Sec. 9), I briefly outline other developments in the theory 
of generalized parton distributions and their applications. 

2 "Old" Parton Distributions 

2.1 Parton Distribution Functions 

The parton distribution function f a /H{ x ) gives the probability that a fast- 
moving hadron H, having the momentum p, contains a parton a carrying the 
momentum xp and any other partons X (spectators) carrying together the 
remaining momentum (1 — x)p. Schematically, 

"X" 



G 
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where the summation is over all possible sets of spectators and ip{H — ► aX} is 
the probability amplitude for the splitting process H — > aX . The summation 




Figure 1: Parton distribution function. 

over X reflects the inclusive nature of the description of the hadron structure 
by the parton distribution functions f a /ii( x )- The parton distribution func- 
tions fa/ii{ x ) have been intensively studied in experiments on hard inclusive 
processes for the last 30 years. The classic process in this respect is the deep 
inelastic scattering (DIS) eN — * e'X whose structure functions are directly 
expressed in terms of f a /ii( x )- 

The standard approach is to write the DIS structure function as the imag- 
inary part of the virtual forward Compton amplitude T^ u (q,p). In the most 
general nonforward case, the virtual Compton scattering amplitude is derived 
from the correlation function of two electromagnetic currents J M (x), J v {y) 

T MV = ij d^x J d 4 ye-^ + ^'y*>(p'\T{J^x)My)}\p) ■ (1) 

In the forward limit, the "final" photon has the momentum q coinciding with 
that of the initial one. The momenta p,p' of the initial and "final" hadrons also 
coincide. The light-cone dominance of the virtual forward Compton amplitude 
is secured by high virtuality of the photons — q 2 = Q 2 > 1 GeV 2 and large total 
center-of-mass (cms) energy of the photon-hadron system s — (p + q) 2 . The 
latter should be above resonance region, with the Bjorken ratio xsj = Q 2 /2{pq) 
fixed. 

An efficient way to study the behavior of Compton amplitudes in the 
Bjorken limit is to use the light-cone expansion for the product 



n MW (x, y) = ¥TJ ll {x)J v {y) 
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of two vector currents in the coordinate representation. The leading order 
contribution is given by two "handbag" diagrams, 

U^(z | X) = -^^s^O^z | X) - ie^ a 5r7 {z | X) J , (2) 

where Sp pva = g^pQvo - gp V g pa + gpag» pi X = (x + y)/2, z = y-x, and 

O a {z\X) = I [$(X - z/2)^{X + z/2) - (*->-*)] , 

<D 5rT (z\X) = 1 [$(X - z/2) lal5 ^(X + z/2) + (z -> -z)] . (3) 

Formally, the parton distribution functions provide parametrization of the 
forward matrix elements of quark and gluon operators on the light cone. For 
example, in the parton helicity averaged case (corresponding to the vector 
operator O a ) the a-quark/antiquark distributions are defined by 

(p | M-z/2)zE(0,z;A)ip a (z/2) | p) 

= u(p)zu{p) J (e-"^f a (x) - e^f- a (x)) dx + 0(z 2 ), (4) 

where E(0, z; A) is the standard path-ordered exponential (Wilson line) of the 
gluon field A which secures gauge invariance of the nonlocal operator. In what 
follows, we will not write it explicitly. Throughout, we use the "Russian hat" 
notation z = "f^z 11 . 

The non-leading (or higher-twist) 0(z 2 ) terms in the above representation 
soften the light-cone singularity of the Compton amplitude, which results in 
suppression by powers of 1/Q 2 (see Sec. 7 for a discussion of twist decompo- 
sition and higher-twist corrections). 

The exponential factors exp[=fix(pz)] accompanying the quark and anti- 
quark distributions reflect the fact that the field ip(z/2) appearing in the op- 
erator ip(—z/2) . . . ip(z/2) consists of the quark annihilation operator (a quark 
with momentum xp comes into this point) and the antiquark creation opera- 
tor (i.e., an antiquark with momentum xp goes out of this point). To get the 
relative signs with which quark and antiquark distributions appear in these def- 
initions, we should take into account that antiquark creation and annihilation 
operators appear in ip(—z/2) . . .ip(z/2) in opposite order. 

In a similar way, one can introduce polarized quark densities A/ a 5 (x) 
which parametrize the forward matrix element of the axial operator 
$ a {-z/2)zj 5 E(0,z;A)ip a (z/2). 

Combining the parametrization (Q) with the Compton amplitude (||) one 
obtains the QCD parton representation 

T»»(p, l)=Y,j Q ^ CteP. l) dx { l + 0(1/Q 2 )| (5) 



8 



Handbook of QCD / Volume 2 




Figure 2: a) General Compton amplitude; 6) s-channel handbag diagram; c) ii-channel 

handbag diagram. 



for the hadron amplitude in terms of the perturbatively calculable hard parton 
amplitude t£ v (xp,q) convoluted with the parton distribution functions f a (x) 
(a = u,d, s, g, . . .) which describe/parametrize nonperturbative information 
about hadronic structure. The short-distance part of the handbag contribution 
is given by the hard quark propagator proportional to l/[(xp + q) 2 + ie]. Its 
imaginary part contains the 5(—Q 2 +2x(qp)) factor (terms 0(p 2 ) are neglected) 
which selects just the x — xsj value from the x-integral. As a result, the DIS 
cross section is directly expressed in terms of f a /Hi x Bj)- Note, however, that 
the factorized representation is valid for the full Compton amplitude in which 
the parton densities are integrated over x. In other words, the variable x in 
Eq.(|) has the meaning of the momentum fraction carried by the parton, but 
it is not equal to the Bjorken parameter XBj- 

The basic parametrization (Q) can be also written as an integral from — 1 
to 1 with a common exponential 

(p | M~z/2)z^ a {z/2) \p)\= u(p)zu(p) J ^ e~ ix ^f a {x) dx + 0(z 2 ). (6) 

For x > 0, the distribution function f a (x) coincides with the quark distribu- 
tion f a (x), while for x < it is given by (minus) the antiquark distribution 
-/«(-*)• 

2.2 Distribution amplitudes 

To give an example of another important type of nonperturbative functions 
describing hadronic structure, namely, the distribution amplitudes, let us con- 
sider the 7*77T° form factor. It is usually measured in e + e~ collisions, but 
for our purposes it is more convenient to represent it through the process in 
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Figure 3: Handbag diagrams and parton picture, a) Virtual forward Compton amplitude 
expressed through the usual parton densities f(x). b) Form factor 7*777° written in terms 
of the distribution amplitude <p(a). 



which electron is scattered with large momentum transfer q off the pion target 
producing a real photon in the final state. The pion, in particular, can belong 
to the cloud surrounding a nucleon. In this case, the j*ir° —> 7 subprocess is 
a part of deeply virtual Compton scattering (DVCS) which will be considered 
in more detail later on. 

The leading order contribution for large Q 2 is given by two handbag dia- 
grams, and in the coordinate representation one deals with the same Compton 
amplitude (||). The only difference is that the nonlocal operators should be 
sandwiched between the one-pion state | ir(r) ) (r is the pion momentum) and 
the vacuum ( |. Since the pion is a pseudoscalar particle, only the axial non- 
local operator 0§ a contributes, and the pion distribution amplitude ip n (a) is 
the function parametrizing its matrix element 

(0 I 4>(-z/2)z l5 ^(z/2) 1 7r(r) > = (rz) J ^ e~ ta(rz ^ n (a)da + 0(z 2 ). (7) 

One can interpret ipn(&) as the probability amplitude 

*{vr(r) -» q(yr) + q((l - y)r)} 

to find the pion in a quark-antiquark state with the pion momentum r shared 
in fractions y = (1 + a)/ 2 and 1 — y = (1 — a)/2. Since the function (p n (a) 
is even in a, the use of the relative fraction a has some advantages when the 
symmetry properties are concerned. 



3 Double distributions 
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3.1 DVCSandDDs 

Now, let us consider deeply virtual Compton scattering (DVCS), an exclu- 
sive process j*(q)N(p) — > "f(q / )N(p') in which a highly virtual initial photon 
produces a real photon in the final state. The initial state of this reaction 
is in the Bjorken kinematics: Q 2 = —q 2 and (pq) are large, while the ratio 
XBj = Q 2 /2{pq) is fixed, just like in DIS. An extra variable is the momentum 
transfer r = p — p 1 . The simplest case is when t = r 2 is small. This does 
not mean, of course, that the components of r must be small: to convert a 
highly virtual photon into a real one, one needs r with a large projection on q. 
Indeed, from q' 2 = it follows that 



i.e. (qr) rts —Q 2 /2 for small t. In the t — > limit, wc can write (qr) = 
xsjiqp)- Taking the initial momentum p in the light cone "plus" direction 
and the momentum q' of the final photon in the light cone "minus" direction, 
we conclude that the momentum transfer r in DVCS kinematics must have a 
non-zero plus component: r + = (p + , with ( = XBj- 

For large Q 2 , the leading order contribution is given again by handbag 
diagrams. A new feature is that the nonperturbative part is described by 
nonforward matrix elements (p — r | . . . \ p) of O a and 0^ a operators. These 
matrix elements are parametrized by generalized parton distributions (GPDs). 
It is instructive to treat GPDs as kinematic "hybrids" of the usual parton 
densities f(x) and distribution amplitudes f(a). Indeed, f(x) corresponds to 
the forward limit r = 0, when the momentum p flows only in the s-channel and 
the outgoing parton carries the momentum xp + . On the other hand, if we take 
p = 0, the momentum r flows in the i-channel only and is shared in fractions 
(l + a)r + /2 and (1 — a)r + /2. In general case, we deal with superposition of two 
momentum fluxes: the plus component k + of the parton momentum k can be 
written as xp + + (1 + a)r + /2. To fully incorporate the symmetry properties 
of nonforward matrix elements, it is convenient to introduce the symmetric 
momentum variable P = (p + p')/2 and write the parton momentum as 



(q + r) 2 = -Q 2 + 2(qr)+t = Q, 



(8) 



k+ = xP+ + (1 + a)r + /2. 



(9) 



This decomposition corresponds to the following parametrization 



(P - r/2|^ a (-z/2)# a (z/2)|P + r/2) 




+0(r) terms + 0(z 2 ), 



(10) 
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xP+(l+a) r/2 



P+r/2 




f(x,a; t) ^P-r/2 



a) 




Figure 4: a) Parton description of deeply virtual Compton scattering in terms of double 
distributions f(x,a;t). b) Support region for f(x,a\t). 

where f a (x,a;t) is the double distribution For the moment, we 

do not write explicitly "O(r)" terms corresponding to u(p')o~ a j3r a z /3 u(p) and 
(rz)u(p')u(p) structures. 



3.2 General properties of DDs 

The support area for f a (x 7 a; t) is shown on Fig. ^b. For any Feynman diagram, 
the spectral constraint |x| + \a\ < 1 can be proved in the a-representationlHI 
using the approach of Ref. 26. Comparing Eq. (|j) with the r = limit 
of the DD definition © gives the "reduction formulas" relating the double 
distribution f a {x, a; t = 0) to the quark and antiquark parton densities 



fa{x,a\t 



>a da = f a (x) 



l+x 



l-\x\ 



i+M 



f a (x, a; t = 0)\ x<a da = -fa(-x) . 



(11) 
(12) 



Hence, the positive- a; and negative- a; components of the double distribution 
f a (x,a;t — 0) can be treated as nonforward generalizations of quark and an- 
tiquark densities, respectively. The usual "forward" densities f a (x) and fa{x) 
are thus given by integrating f a (x,a; t — 0) over vertical lines x — const for 
x > and x < 0, respectively. In principle, we cannot exclude the third possi- 
bility that the functions f a {x, a; t) have singular terms at x = proportional to 
5(x) or its derivative (s). Such terms have no projection onto the usual parton 
densities. We will denote them by fM(x,a;t) — they may be interpreted as 
coming from the ^-channel meson-exchange type contributions. In this case, 
the partons just share the plus component of the momentum transfer r: in- 
formation about the magnitude of the initial hadron momentum is lost if the 
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exchanged particle can be described by a pole propagator l/(t — m 2 M ). Hence, 
the meson-exchange contributions to a double distribution may look like 

/+(x,y;t)~5(x)-4^: or fa (x, a; t) ~ 5'{x) ^fM. , e tc. , (13) 

m M m M 

where ip^i 01 ) are the functions related to the distribution amplitudes of the 
relevant mesons M ± . The two examples above correspond to x-even and £-odd 
parts of the double distribution f(x,a;t). 

Due to hermiticity and time-reversal invariance properties of nonforward 
matrix elements, the DDs are even functions of a, 

f a {x,a;t) = f a (x,-a;t) . 

In particular, the functions <£>^-(a:) for singular contributions f^(x,a;t) arc 
even functions ip^ (a) = tp^ (— a) of a both for x-even and x-odd parts. 

Note that the \i <-> v symmetric part of the DVCS amplitude contains only 
the C-even operators O a a . Their matrix elements are parametrized 

(p',s'\ z°O a (-z/2,z/2) \p,s)\ z 2 =0 

r l f l-\x\ 

= u{p',s')zu(p,s) dx e- ix ( Pz ^- ia ( rz V 2 f%(x,a; t)da 

J-l J-l+\x\ 

+0(r) terms (14) 
by the DDs 

fa{x,<x; i)=sign(a;)[/a(|a;|,a; t) + fa{\x\, a ; t)] 

which are odd functions of a;. In applications to the hard meson electropro- 
duction one also needs valence-type DDs 

fa(x,a; t) = [/ a (|o;|,a; t) - f a (\x\,a; t)] 

parametrizing matrix elements of C-odd operators ■ip a (—z/2)ztp a (z/2) + 
$ a {z/2)z*p a {-z/2). 



4 Skewed parton distributions 

4-1 General definition 

An important parameter for nonforward matrix elements is the coefficient of 
proportionality Q = r + /p + (or £ = r + / P + ) between the plus components of the 
momentum transfer and the initial (or average) hadron momentum. It specifies 
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the skewedness of the matrix elements. The two skewedness parameters are 
related by 

*=23C" < 15 > 
The characteristic feature implied by the definition of the double distribu- 



tion (10) is the absence of the ^-dependence in f(x, a;t). An alternative way 
to parametrize nonforward matrix elements of light-cone operators is to use £ 
and the total momentum fraction x = x + £ a as independent variables and 
introduce skewed parton distributions (SPDs). The shape of SPDs explicitly 
depends on the skewedness of the relevant nonforward matrix element. 

There are two types of SPDs: off-forward parton_distributions0Q (OF- 
PDs) H(x,£;t) and nonforward parton distributionsBo (NFPDs) T$(X',t). 
The basic difference is that the skewedness parameter £ (£) and the parton 
momentum k + in the OFPD (NFPD) formalism is measured in units of the 
average (initial) momenturtuP + (p + ). Hence, there are one-to-one relations be- 
tween OFPDs and NFPDstZl We start with OFPDs because they have simpler 
symmetry properties. 

The relation between OFPDs H(x , £; t) and DDs f(x, a; t) is given just by 
the change of variables from {x, a} to {x, £}, 

,1 .l-|x| 

H(x,£;t) = / dx 5(x + £a — x) f(x, a; t) da . (16) 

J-l J-l+\x\ 

If we require that the light-cone plus components of both the momentum 
transfer r and the final hadron momentum p — r are positive (which is the 
case for DVCS), then < £ < 1 and < £ < 1. Using the spectral property 
< |x| + |a| < 1 of double distributions we obtain that the OFPD variable x 
satisfies the constraint < \x\ < 1. Note also that Eq. ([l6]) formally allows to 
construct H(x,^; t) both for positive and negative values of £. Since the DDs 
f(x,a;t) are even functions of a, the OFPDs H(x, are even functions of 

H(x,t;t)=H(x,-Z;t). 

This result was originally obtained by X. Jill! who used hermiticity and time 
reversal invariance properties in the direct definition of OFPDs 

{p'\M-z/2)zip a (z/2)\p) 

= u(p')zu(p)J e^ Pz ^ H a (x^;t)dx + 0{r) term +0(z 2 ). (17) 
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xP+(l+a)r/2 xP-(l-a)r/2 a 




a) b) c) 



Figure 5: Description of the nonforward matrix element in terms of a) double distribution 
f{x,a\t) and b) off-forward parton distribution H(x,£;t). c) Integration lines for integrals 

relating OFPDs and DDs. 

4-2 Structure of SPDs 

The parton interpretation of the OFPD definition is the following: the quark 
going out of the parent hadron carries the fraction x + £ of the average hadron 
momentum P + , while the momentum of the "returning" quark is (pt — £)-P + . 
For definiteness, we shall assume below that £ is positive. Treating a quark 
with a negative momentum as an antiquark, we can distinguish 3 components 
of H(x, £;t). For x > £ both the initial fraction x + £ and the final one 
x — £ arc positive. Hence, H(x, in this region corresponds to a modified 
quark distribution. Similarly, for x < — £ both fractions are negative, and 
H(x^]t) can be treated as an antiquark distribution. In the third (middle) 
region, — £ < x < £, the two fractions have opposite signs, and i/(a;,£;£) 
describes splitting of a quark-antiquark pair from the initial hadron. The total 
momentum carried by the qq pair is 2£P + = r + , and it is shared in fractions 
(x + 0P + = (1 + P)r + /2 and (£ - x)P+ = (1 - 0)r+/2, where = x/£. 
Note that \(3\ < 1 in the \x\ < £, region. Hence, the third component can be 
interpreted as the probability amplitude for the initial hadron with momentum 
(1 + £)P + to split into the final hadron with momentum (1 — and a two- 

parton state with total momentum r + shared by the partons of the pair. Thus, 
we may expect that H (x,£; t) in the middle region looks more like a distribution 
amplitude. 

The relation between DDs and SPDs can be illustrated on the DD support 
rhombus |x| + |a| < 1 (see Fig. ^p). The delta-function in Eq. ( p^ ) specifies the 
line of integration in the {x, a} plane. To get H(x,^;t), one should integrate 
f{x, a) over a along a straight line x — x — £a. Fixing some value of £, one 
deals with a set of parallel lines intersecting the x-axis at x — x. The upper 
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limit of the a-integration is determined by intersection of this line either with 
the line x + a = 1 (this happens if x > £) or with the line —x + a = 1 (if x < £). 
Similarly, the lower limit of the a-integration is set by the line x — a = 1 for 
x > — £ or by the line x + a = —1 for x < — £. The lines corresponding to 
5 = ±£ separate the rhombus into three parts generating the three components 
of H(x,£;t): 

H a (x,£;t)=0(£ <x< 1) / 1 ' f a (x-£a,a)da 
+0(-£ <x<Q f a (x~ fa, a) rfa 

/ 1 — S 

— "TTF 

l + £ 

+0(-l < £ < -0 / 1+ ' / a (x - fa, a) da. (18) 

Recall that integrating the DD /(a;, a; i = 0) over a vertical line gives the 
usual parton density /(&). To get the t = SPDs one should scan the same 
DD along the lines having a ^-dependent slope. Thus, one can try to build 
models for SPDs using information about usual parton densities. Note, how- 
ever, that the usual parton densities are insensitive to the meson-exchange type 
contributions Hm(x, £; t) coming from the singular x — parts of DDs. Thus, 
information contained in SPDs originates from two physically different sources: 
meson-exchange type contributions iJ/vf(i,f;i) and the functions Hm(x, f;i) 
obtained by scanning the x ^ parts of DDs f(x,a;t). The support of ex- 
change contributions is restricted to |x| < f. Up to rescaling, the function 
#Af(£,f;i) has the same shape for all f, e.g., <pm(x/£; £)/|£|. For any nonva- 
nishing f , these exchange terms become invisible in the forward limit f — ► 0. 
On the other hand, interplay between x and f dependences of the component 
resulting from integrating the x ^ part of DDs is quite nontrivial. Its support 
in general covers the whole — 1 < x < 1 region for all f including the forward 
limit f in which they convert into the usual (forward) densities f a (x), f a {x). 
The latter are rather well known from inclusive measurements, at small t. 



4-3 Polynomiality and analyticity 

In our derivation, DDs are the starting point, while SPDs are derived from 
them by integration. However, even if one starts directly with SPDs, the 
latter possess a property which can be incorporated only within the formalism 
of double distributions. Namely, the x N moment of H(x,£ i ;t) must be an 
Nth. order polynomial of f . This restriction on the interplay between x and 
f dependences of H(x,£;t) followstj from the simple fact that the Lorentz 
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indices fix . . . /ijy of the nonforward matrix element of a local operator O^i - pn 
can be carried cither by P Mi or by r Mi . As a result, 



(P - r/2\m(d + ) N m\P + r/2) = E U ) (P + ) N ~ k {r+) k A 



k=0 



y* e; 5 w dx = f2 c k y * ^ y ' 



N / \ 
fc=0 ^ ' 

where (^) = N\/(N — k)\k\ is the combinatorial coefficient. Our derivation 
(16) of OFPDs from DDs automatically satisfies the polynomiality condition 
(19), since 

N /at\ rl 

f{x,a)x N - k a k da . (20) 

\ /. II: I ill 

k=0 

Hence, the coefficients Ajsk in (|l^) are given by double moments of DDs. This 
means that modeling SPDs one cannot choose the coefficients A^k arbitrarily: 
symmetry and support properties of DDs dictate a nontrivial interplay between 
N and k dependences of An^s. After this observation, the use of DDs is a 
necessary step in building consistent parametrizations of SPDs. 

The formalism of DDs also allows one to easily establish some important 
properties of skewed distributions. Notice that due to the cusp at the upper 
corner of the DD support rhombus, the length of the integration line nonanalyt- 
ically depends on x for x = ±£. Hence, unless a double distribution identically 
vanishes in a finite region around the upper corner of the DD support rhombus, 
the x-dependence of the relevant nonforward distribution must be nonanalytic 
at the border points x — ±£. Still, the length of the integration line is a con- 
tinuous function of x. As a result, if the double distribution f(x,a;t) is not 
too singular for small x, the skewed distribution H (x, £; t) is continuous at the 
nonanalyticity points x = ±£. Because of the l/(i ± £) factors contained in 
hard amplitudes, this property is crucial for PQCD factorization in DVCS and 
other hard electroproduction processes. 

Note, that there may be also the exchange contributions Hm(%, £> t). If it 
comes from a 8(x)ip(a) type term and ip(a) vanishes at the end-points a = ±1, 
the ffjtf (a;,£;t) part of SPD vanishes at x = ±£. The total function H(x, £;t) 
is then continuous at the nonanalyticity points x = ±£. In the C-even case, 
the DDs should be odd in x, hence the singular term involves S'(x)ip(a) (or 
even higher odd derivatives of 5(x)). One can get a continuous SPD in this 
case only if tp'(a) vanishes at the end points. 
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4-4 Nonforward parton distribution functions 

In the NFPD formalism!*! the skewedness parameter £ = r + /p + and the parton 
momentum k + are measured in units of the initial momentum p + . Again, 
one can start with double distributions F(x,y;t) writing the outgoing parton 
momentum as k + = xp + + yr + and that of the returning one as xp + — (1 — y)r + 
(see Fig. [|a). The support area for F(x,y;t) is specified by — 1 < x < 1, 
< y < 1, < x+y < 1 (see Fig. ||c). The relation between F a (x, y; t = 0) and 
the usual quark and antiquark densities is given by the "reduction formulas" 

l-X j-l 

F a (x, y; t = 0)\ x>0 dy = f a (x) ; / F a (x, y; t = 0)| x<0 dy = -/ a (-x) . 

J —x 

(21) 

If we define the "tilded" DDs by 

F a (x, y; t) = F a (x, y; t) x>0 ; F a (x, y; t) = -F a (-x, 1 - y; t) x<a , 
then x is always positive and the reduction formulas 



Fa : a {x, y,t 



dy = fa,a( X ) 



(22) 



have the same form in both cases. The new antiquark DDs Fa(x 7 y]t) also 
"live" on the triangle < x, y, x + y < 1. 



xp+yr 



I \ xp-(l-y) r 



(x-Vp 





C X 



c) 



Figure 6: Description of nonforward matrix elements in terms of a) double distributions 
F(x,y;t) and b) nonforward parton distributions JF*(X;t). c) Integration lines for the inte- 
grals relating NFPDs and DDs. 

The a — > —a symmetry of the OFPD-oriented DDs f(x,a;t) corresponds 
to the symmetry of F a ^(x, y; t) with respect to the interchange y <-> 1 — x — y 
("Munich symmetry") established in Ref. 27. 

Using 7' + = (p + and introducing the total momentum fraction X = x + (y, 
we define the nonforward parton distributions (see Fig. ^b) 

r c (X;t)= [ dx f 5(x + (y-X)F t (x,y;t)dy . (23) 
Jo Jo 
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The relation between NFPDs and DDs Fi(x,y) can be illustrated on the 
DD support triangle < x,y,x + y < 1 (see Fig. |6p). To get !F^(X;t), one 
should integrate F(x, y; t) over y along a straight line x — X — C,y. The upper 
limit of the y-integration is determined by intersection of this line either with 
the line x + y = 1 (this happens if X > £) or with the y-axis (if X < £) : 

f (i-x)/(i-Q 

FAX-t) = 8(X>() F t (X-(y iy ;t)dy 



Jo 

r-x/C 

+ 9(X<0 / F l (X-(y,y;t)dy. (24) 
Jo 

The returning parton carries the fraction X — £ of the initial hadron mo- 
mentum p + , which is positive in the X > £ region (where NFPDs can be 
treated as modified parton densities) and negative in the X < £ region (where 
NFPDs resemble distribution amplitudes). 

For reference purposes, we present the relations between the NFPD and 
OFPD variables 

C -TT~C ~ IT?' X=X ~ C ~ITI' x ~T^j2- (25) 

4-5 Relation to form factors 

GPDs depend on the invariant momentum transfer t, hence, they may also be 
treated as generalizations of hadronic form factors. In particular, H a (x, £; t) is 
related to the Dirac form factor of the proton!^ 



e a J H a (x,£,t)dx = F 1 (t), (2(i) 



where e a is the electric charge of the relevant quark. Just like for form factors, 
there are extra generalized parton distributional^ E a (x,£,t) corresponding to 
hclicity flip in the nonforward matrix element. They are related to F 2 (t) form 
factor, 

>„ / E a (x,Z,t)dx = F 2 (t). (27) 



E- 



Note that though the shape of GPDs changes when £ is varied, the integrals 
(p6|), ( p7| ) do not depend on £. Since the E functions are accompanied by the 
r M factor in the parametrization of the nonforward matrix element 

( P '\i>j^\p) ~ u(p') ltl u( P yH» + J^u{j>')a^r v u{pYE^ , 
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they arc invisible in deep inelastic scattering described by exactly forward r = 
Compton amplitude. However, the t = 0,£ = limit of the E distributions 
exists, 

E a (x,Z = 0;t = 0)=e a (x) . 

In particular, for t = the integral (|2^) gives the anomalous magnetic moment. 
Moreover, the recent interest to DVCS aniL generalized parton distributions is 
largely due to observation made by X. JiEl that the integral 

1 f 1 

J i = oYl / \J«( x ) + e ^ xdx ( 28 ) 

is related to the total spin and orbital momentum contribution of the quarks 
into the proton spin. Ji proposed to use deeply virtual Compton scattering to 
get access to the e a (x) functions. These function can be also accessed in hard 
meson electroproduction processesE3 

The DVCS amplitude contains two other generalized parton distributions 
H(x,£,;t) and E(x- r £]t). They parametrize the nonforward matrix element of 
the axial operators 

(p'l^vwflp) ~ uip'hMp) + ^(p'Hp) ■ (29) 

In the forward limit, the iJ-distributions reduce to the polarized parton 
densities H a (x,£ = 0;t = 0) = Af a (x). After the x-integration, the 
H a (x,t; = 0;t = 0) distributions produce the flavor components of the axial 
form factor GU(0- Similarly, the functions E(x,^;t) are related to the pseu- 
doscalar form factor Gp(t). At small t, they are dominated by the pion pole 
term l/(t— m%). 

4-6 Gluon distributions 

In a similar way, we can write parametrizations defining double and skewed 
distributions for gluonic operators 

(P - r/2 | z^ Zu G; a (-z/2)E ab (-z/2, z/2; A)G h av {zl2) \ P + r/2) 

r l f l-\x\ 

= (Pz)u(p')zu{p) dx e-^ Pz ^ ta< - rz) / 2 xf g (x,a;t)da + ... 

J-l J-l+\x\ 

= (Pz)u(p')zu(p) J e~ lS:{Pz) H g (x^-t)dx + ... . (30) 

Note, that the gluon SPD H g (x, £; t) is constructed from xf g (x, a; t). Just like 
the singlet quark distribution, the gluon double distribution f g (x,a;t) is an 
odd function of x. 
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5 Modeling GPDs 

There are two approaches used to model GPDs. One is based on a direct calciu 
lation of parton distributions in specific dynamical models, such as baftpjadelpj 
chiral soliton modelpl light-cone formalismpl etc. Another approacli^OEl is a 
phcnomenological construction based on reduction formulas relating GPDs to 
usual parton densities f(x), A/(x) and form factors Fi(t), i*2(i), GU(i), Gp(t). 
The most convenient way to construct such models is to start with double 
distributions f(x,a;t). 



5.1 Modeling DDs 

Let us consider the limiting case t = 0. Our interpretation of the x- variable as 
the fraction of the P momentum and the reduction formula stating that the 
integral of f a (x,a) over a gives the usual parton density f a (x) suggest that 
the profile of f a (x,a) in the x-direction follows the shape of f a {x). Thus, it 
make sense to write 

f(x,a) = h(x,a)f(x), (31) 
where the function h(x, a) normalized by 

1 — |a?| 

h(x,a)da = 1 (32) 

-1+1*1 

characterizes the profile of f(x, a) in the a-direction. The profile function 
should be symmetric with respect to a — > ~a because DDs f(x, a) are even in 
a. For a fixed x, the function h(x, a) describes how the longitudinal momentum 
transfer r + is shared between the two partons. Hence, the shape of h(x, a) 
should look like a symmetric meson distribution amplitude <p(a). Recalling 
that DDs have the support restricted by \a\ < 1 — |x|, to get a more complete 
analogy with DAs, it makes sense to rescale a as a = (1 — \x\)/3 introducing the 
variable (3 with x-independent limits: — 1 < f3 < 1. The simplest model is to 
assume that the /3-profile is a universal function g(0) for all x. Possible simple 
choices for g((3) may be S(/3) (no spread in /3-direction) , /3 2 ) (characteristic 
shape for asymptotic limit of nonsinglet quark distribution amplitudes), y§(1 — 
/3 2 ) 2 (asymptotic shape of gluon distribution amplitudes), etc. In the variables 
x, a, this gives 

h^(x,a)=S(a) , h<H*,«) = l {1 ~Wy-/ , 
hWt \ 15 [(1 - \x\) 2 - a 2 } 2 

fe (x ' a) ^i6 (i-NP ■ (33) 
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These models can be treated as specific cases of the general profile function 
h*Hx a) = T(2b + 2) Kj-N)!^3! (34) 

2 2b +!r 2 (6 + i) (i-|z|) 2b+1 ' { ' 

whose width is governed by the parameter 6. 
5.2 Modeling SPDs 

Let us analyze the structure of SPDs obtained from these simple mod- 
els. In particular, taking f(°°\x,a) = 5(a) f(x) gives the simplest model 
ij(°°J (J, £; t = 0) = f(x) in which OFPDs at t = have no ^-dependence. For 
DDs producing nonforward parton distributions !F^(X;t — 0), this is equiva- 
lent to the F(°°\x, y) = 5(y — x~/2)f(x) model, which gives 

i.e., NFPDs for non-zero C are obtained from the forward distribution f(X) = 
Tq = o(X) by shift and rescaling. 

In case of the 6=1 and 6 = 2 models, simple analytic results can be 
obtained only for some explicit forms of f(x). For the "valence quark" -oriented 
ansatz f^(x,a), the following choice of a normalized distribution 

/ (1) W = ^8p\^°(l-^) 3 (36) 
6 1 (1 — a) 

is both close to phenomenological quark distributions and produces a simple 
expression for the double distribution since the denominator (1 — x) 3 factor in 
Eq. ( |33| ) is canceled. As a result, the integral in Eq. ( |l8| ) is easily performed 
and we get 

H^(x,0\ m >s = ^(l-^){[(2-aK(l-i)(4" a + ^ a ) 

+ (£ 2 - x)(x 2 - a - x 2 _- a )] 9(x) + (x^ -£)} (37) 



for \x\ > £ and 

HU^ 0I|*|<£ = ^ (l - |) R~ Q [(2 - <0£(1 - x) + (e - x)] + (x^ -x)} 

(38) 

in the middle — £ < x < £ region. We use here the notation x + = 

and X- — (x — £)/(l — £)■ To extend these expressions onto negative values 

of £, one should substitute £ by |£|. One can check, however, that no odd 
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powers of |£| would appear in the x moments of H 1V (i,£). Furthermore, 
these expressions are explicitly non-analytic for x — ±£. This is true even if 
a is integer. Discontinuity at x = ±£, however, appears only in the second 
derivative of H iv (x,£,), i.e., the model curves for H iv {x : t;) look very smooth 
(see Fig. where the curves for NFPDs are also shown) . 




0.2 0.4 0.6 0.8 1 -1 -0.5 . 5 1 



Figure 7: Valence quark distributions: NFPDs F?(X) (left) and OFPDs H^(x,£) (right) 
with a = 0.5 for several values of f, namely, 0.1, 0.2, 0.4, 0.6, 0.8 and corresponding values 
of ^ = C/(2 — £)■ Lower curves correspond to larger values of f. 



For a = 0, the x > £ part of OFPD has the same x-dependence as its 
forward limit, differing from it by an overall ^-dependent factor only, 

H 1V (x, OU=o = 4 > + 2 9(\x\ < . (39) 

The (1 — |i|) 3 behavior can be trivially continued into the \x\ < £ region. How- 
ever, the actual behavior of H iv (x, £)| Q= o hi this region is given by a different 
function. In other words, H 1V (x, £)| a =o can be represented as a sum of a func- 
tion analytic at border points and a contribution whose support is restricted 
by |*| < It should be emphasized that despite its DA-like appearance, this 
contribution should not be treated as an exchange-type term. It is generated 
by the regular x ^ part of the DD, and, unlike the <p(i/£)/£ functions its 
shape changes with £, the function becoming very small for small £. 

For the singlet quark distribution, the DDs f s (x,a) should be odd func- 
tions of x. Still, we can use the model like (^6|) for the x > part, but take 
f s (x, a)\ x ^o = A f^(\x\, a) sign(x). Note, that the integral ( |l8| ) producing 
H s {x,Q in the \x\ < £ region would diverge for a — > x/£ if a > 1, which 
is the usual case for standard parametrizations of singlet quark distributions 
for sufficiently large Q 2 . However, due to the antisymmetry of f s (x,a) with 
respect to x —> —x and its symmetry with respect to a — > —a, the singularity 
at a = x/£ can be integrated using the principal value prescription which in 
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-0.75 - 




Figure 8: Model for singlet quark distributions: OFPDs ffi(5;,£) for values of f corre- 
sponding to £ equal to 0.2, 0.4, 0.6. Lower curves correspond to larger values of 

this case produces the x — > — a; antisymmetric version of Eqs. (|37]) and (pq). 
For a = 0, its middle part reduces to 

0||3|<«,«=o = 2* 3 ^7 ( i+^ a;2 • ( 4 °) 
The shape of singlet SPDs in this model is shown in Fig. || 

5.3 Poly akov- Weiss terms 

Note, that the operator -ip(—z/2)zip(z/2) is proportional to z. Hence, 
parametrizing its matrix elements in terms of parton distributions, it makes 
sense to use the structures which are also linear in z, like u(p')zu(p), 
u(p')(zf — rz)u(p) for the nucleon target, (Pz) for the pion target, etc. How- 
ever, there is a subtlety emphasized by Polyakov and WeissEl Namely, using 
parametrization by DDs, one treats (Pz) and (rz) as independent variables. 
This means that in the pion case, e.g., one should deal both with (Pz) and 
(rz) structures: 

(P - r/2 | ^(-z/2)zil){z/2) \ P + r/2) 
= 2(Pz) f dx f '' e -**(P*)-to(r*)/a f(x,a;t) da 

J -l J—l+\x\ 

+ (rz) J e- ia( - rz ^ 2 D(a;t)da , (41) 

where the Polyakov- Weiss (PW)§ term D(a;t) accumulates the r Ml . . . r^ n 
parts of the expansion for the matrix elements of the local operators 
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i/j(0)7 Ml D M2 . . . D^ n ip(0). There is no sensitivity in the D-term contribution to 
the value of the average momentum P term: the parton momenta depend only 
on r. Hence, the .D-term is a particular case of the exchange contributions. 

Switching to skewed parton distributions, one deals with just one structure 
(Pz), and one can incorporate the PW term contribution!!^ D(x/0, t)9(\x\ < 
into H(x, £;t). 

In the nucleon caseP the additional structure is (rz)u(p')u(p). As a result, 
the skewed distributions H, E have two components, one is obtained from the 
relevant DD and another comes from the D-term 

H(x, e t) = H DD (x, C; t) + D{x/0, t)B(\x\ < , (42) 
E(x,0,t) = E DD (x,0,t) - D(x/0,t)6(\x\ < . (43) 

Note, that the D-term drops from the Ji's sum rulc@@ since H + E = Hd d + 
Edd- 

It should be noted that explicit calculations of skewed parton distribu- 
tions performed within the chiral soliton modeO show that the middle region 
behavior of SPDs strongly resembles that of distribution amplitudes. 

5.4 Inequalities 

In the case when x > £, the integration line producing H(x, £;t) (see Fig. ||c) 
is inside the space between two vertical lines giving the usual parton densities 
/(x_) and f(x + ), with x+ = (5 + 0/(1+0 and a;_ = (x - 0/(1 - £): 

r x+ dx 

H a {x,Zit) = 0(Z<x<l) f a (x,(x-x)/Z)—dx + ... . (44) 

J X- s 

The combinations X—, x+ have a very simple interpretation: they measure the 
momentum of the initial or final parton in units of the momentum of the rele- 
vant hadron. Assuming a monotonic decrease of the double distribution f(x, a) 
in the ^-direction and a universal profile in the a-direction, one may expect 
that H a (x,£;t) is larger than f(x+) but smaller than /(#_). Inequalities be- 
tween forward and nonforward distributions were discussed in Refs. 31, 19, 32. 
They are based on the application of the Cauchy-Schwartz inequality 

|(jv[(i - on (* - OP; s\ N[(i + op]; (x + OP; s)\ 2 
< £<;v[(i + OP}; {x + OP; s\ N[{i + OP}; (£ + 0^; S) 

s 

x £<jv[(i - OP]; (£ - OP; s'\ n[(i - OP]; {* - OP; s') (45) 

S' 
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to the skewed parton distributions written generically as 

H(x, & t) = ^2(N[(i - OP}; (i - OP; s\ N[(l + op}; (i + OP; S) , (46) 
s 

where \N[(1 ± OP}; i OP; describes the probability amplitude that the 
nucleon with momentum (1 ± £)P converts into a parton with momentum 
(x ±£)P and spectators S. The forward matrix elements give the usual parton 
densities 

£(JV[(i ± OP}; (5 ± 0^; ^) = TT7 /(a;±) - (47) 

As a result, one obtainsSEHEI for the quark distributions 

H a {x,0t) < 
For the gluon distribution, one has 



fa{x+)fa{X-) 



(48) 



Hg(x,0,t) < \J x +X-fg{x+)f g (X-) . 



(49) 



It is clear that the whole consideration makes sense only if x > If 
x < the integration line x — x — t;a intersects the line x = 0, where the 
usual parton densities are infinite. Furthermore, when negative x are involved, 
the behavior of DDs along the line cannot be monotonic. Another deficiency 
of the Cauchy-Schwartz-type inequalities is that they do not give the lower 
bound for nonforward distributions though our graphical interpretation sug- 
gests that H(x, 1 = 0) is larger than f{x+) if the s-dependence of the double 
distribution f(x, a) along the lines a — k(l — \x\) is monotonic. 

5.5 SPDs at small skewedness 

To study the deviation of skewed distributions from their forward counterparts 
for small £ (or £), let us consider the x > £ part of H(x,0 (see Eq. (|l^)) and 
use its expansion in powers of £ 



H(x;0\x>i 



1 f^-V d 2 f(x,a) 

2 J dx 2 

.df(x,a) 



a da 



(1-2) 

df{x,a) _ 
da 



di 



(50) 



where f(x) is the forward distribution. For small £, the corrections are formally 
O(0)- However, if f(x,a) has a singular behavior like x~ a , then 

d 2 f(x,a) a(l + a) 



dx 2 



-f{x,a) , 
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and the relative suppression of the first correction is 0(£ 2 /x 2 ). Though the 
corrections are tiny for i 3> £, in the region x ~ £ they have no parameter 
smallness. It is easy to write explicitly all the terms which are not suppressed 
in the x ~ £ — > limit 

(51) 

where the ellipses denote the terms vanishing in this limit. This result can be 
directly obtained from Eq. ( |l8| ) by noting that for small x, we can neglect the 
x-dependence in the limits ±(1 — |a;|) of the a-integration. Furthermore, for 
small x one can also neglect the x-dependence of the profile function h(x, a) in 
Eq. ( |3"l| ) and take the model fix, a) = f(x)p(a) with p(a) being a symmetric 
normalized weight function on — 1 < a < 1. Hence, in the region where both 
x and £ are small, we can approximate Eq. ([l8]) by 

H(x; = "P" J ^ f{x - £a)p{a) da + . . . , (52) 

i.e., the OFPD H(x; £) is obtained in this case by averaging the usual (forward) 
parton density fix) over the region i — £<x<x + £ with the weight p(a). 
The principal value prescription "P" is only necessary in the case of singular 
quark singlet distributions which are odd in x. In terms of NFPDs, the relation 
is 

j£- c (X) = "P"y fiX -(il+a)/2)pia)da + ... , (53) 

i.e., the average is taken over the region X — ( < x < X. 

In fact, for small values of the skewedness parameters £, £, one can use 
Eqs. (|5^), d53| ) for all values of x and X: if x ^> £, Eq. ( ]5^ ) gives the correct 
result -ff(i;£) = /(i) + 0(£ 2 ). Hence, to get SPDs at small skewedness, one 
only needs to know the shape of the normalized profile function p(a). 

The imaginary part of hard exclusive meson electroproduction amplitude is 
determined by the skewed distributions at the border point x = £ (or X = Q. 
For this reason, the magnitude of !F^iQ [or and its relation to the 

forward densities fix) has a practical interest. This example also gives a 
possibility to study the sensitivity of the results to the choice of the profile 
function. Assuming the infinitely narrow weight p(a) = (5(a), we have J-^iX) — 
fiX - C/2) + ... and H(x, £) = fix). Hence, both JF C (£) and , £) are given 
by fixBj/2) because C = XBj and £ = XBj/2 + ■ ■ ■■ Since the argument of fix) 
is twice smaller than in deep inelastic scattering, this results in an enhancement 
factor. In particular, if fix) ~ x~ a for small x, the ratio lZi() = ^(C)//(C) 
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is 2 a . The use of a wider profile function p(a) produces further enhancement. 
For example, taking the normalized profile function 

" fc(a) 55 r(i/2)r(6 + i) (1 - a > = 2^W + D (1 " } (54) 

and /(x) ~ x _a we get 

nW (c) = ^(0 _ r(2fr + 2)r(6-q + l) 

* (C) - /(C) ~r(2&-a + 2)r(& + i) (55) 

which is larger than 2° for any finite b and < a < 2. The 2 a enhancement 
appears as the n — > oo limit of Eq. (pi|). For small integer n, Eq. ( |54| ) reduces 
to simple formulas obtained in Refs. 21, 22. For n = 1, we have 

/(C) (l-a/2)(l-a/3)' 1 J 

which gives the factor of 3 for the enhancement if a = 1. For 6 = 2, the ratio 
becomes 



•7f =2) (0 1 



/(C) (l-a/3)(l-o/4)(l-a/5) 



(57) 



producing a smaller enhancement factor 5/2 for a = 1. Calculating the en- 
hancement factors, one should remember that the gluon SPD ^(X) reduces 
to Xfg(X) in the C = limit. Hence, to get the enhancement factor corre- 
sponding to the f g (x) ~ x _A small-x behavior of the forward gluon density, 
one should take a = A — 1 in Eq. (|54]). As a result, the 1/x behavior of the 
singlet quark distribution gives the factor of 3 for the IZ^ (£) ratio, but the 
same shape of the gluon distribution results in no enhancement. 

Due to evolution, the effective parameter a characterizing the small-x be- 
havior of the forward distribution is an increasing function of Q 2 . Hence, for 
fixed b, the IZ^^) ratio increases with Q 2 . In general, the profile of f(x,a) 
in the a-direction is also affected by the PQCD evolution. In particular, in 
Ref. 21 it was shown that if one takes the ansatz corresponding to an extremely 
asymmetric profile function p{a) ~ 8(1 + a), the shift of the profile function 
to a more symmetric shape is clearly visible in the evolution of the relevant 
SPD. In the next sections, we will discuss the evolution of GPDs and study 
the interplay between evolution of x and a profiles of DDs. 
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6. 1 Evolution kernels for double distributions 

The QCD perturbative expansion for the matrix element in Eq. ( |Tc| ) generates 
at one loop level the terms proportional to lnz 2 . In other words, the limit 
z 2 — > is singular and the distributions F(x,y;t), f(x,a;t), etc., contain 
logarithmic ultraviolet divergences which require an additional i?-operation 
characterized by some subtraction scale /i: F(x,y;t) — ► F(x,y;t \ fi). The 
/i-dependence of F(x, y;t \ fx) is governed by the evolution equation 

fx^-F a {x,y;t\fi) = [ f y2R ab (x,y;C,r l )F b ^,i r ,t\ f i)9^ + f 1 <l)d^dr ! , 
d H Jo Jo V 

(58) 

where a,b — g, q. A similar set of equations, with the kernels denoted by 
AR ab (x,y;£,r]) governs the evolution of the parton helicity sensitive distribu- 
tions G a (x, y;t\{i). Since the evolution kernels do not depend on t, from now 
on we will drop the i-variable from the arguments of F(x, y; 1 1 fx) in all cases 
when this dependence is inessential (likewise, the /i-variablc will be ignored in 
our notation when it is not important). 

Since integration over y converts F a (x,y;t = Q | >u.) into the parton dis- 
tribution function f a (x\fx), whose evolution is described by the Dokshitzer- 
Gribov-Lipatov-Altarelli-Parisi (DGLAP) 33-35 equations 

fx-ff a (x \n)= [ fp ab (x/Oh(t; \ (59) 
dfj, J x £ 

the kernels R ab (x, y; £, rj) must satisfy the reduction relation 

R ab {x,y^,ri)dy = -P ab {x/0. (60) 

Alternatively, integration over x converts F a (x, y;t = | fx) into an object sim- 
ilar to a meson distribution amplitude (DA), so one may expect that the result 
of integration of R ab (x, y; £, rj) over x should be related to the kernels govern- 
ing the evolution of distribution amplitudes [Efremov-Radyushkin-Brodsky- 
Lepage (ERBL) 5 - 6 ' 37 evolution], e.g., in case of the qq kernel 

" K>*{x,y;Z,r))dx = V**{y,Ti). (61) 







These reductionjiroperties of the R qq (x, y; £, rf) kernel can be illustrated using 
its explicit formp] 

B^{x,y;C,mg) = —cMe(0 < x/t < min{y/ V ,y/f}}) 

IT £ I 
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0(0 < x/tj < l)x/j 

(i-^/O 

-6{l-x/0S{y-v) 



Jl 
1 

2 + 



vhi) + ~s(x/£- 

V 



■ dz 



y/v) 



i 



(62) 



Here the last (formally divergent) term, as usual, provides the regularization 
for the l/(x — t;) singularities present in the kernel. This singularity can be also 
written as l/(r) — y) for the term containing 8(x / '£ — y j 'rf) and as l/(fj — y) for 
the term with 8(x/t; — y/ff). Depending on the chosen form of the singularity, 
incorporating the 1/(1 — z) term into a plus-type distribution, one should treat 
z as y/rj or y/fj. One can check that integrating R qq (x,y; £,??) over y or 
x gives the DGLAP splitting function 

'1 + z 2 " 



P qq {z-g) 



C F 



1 - z 



(63) 



and the Brodsky-Lepage evolution kernel 



-C f 



1 



y-r) 



rj-y 
9(y > 77) 



0(y < v) 



(64) 



Here, "+" denotes the standard "plus" regularizationt 



6.2 Light-ray evolution kernels 

The interrelation between different types of evolution kernels follows from the 
fact that, in the leading logarithm approximation, the evolution equations can 
be written for the light-cone operators themselves, 37-40 without any reference 
to particular matrix elements 

7 /• 1 /• 1 

u—OJO,z)= / B ab (u,v)O b (uz,vz)0(u + v < l)dudv, (65) 
Jo Jo ^ 

where v = 1 — v and a,b = q, g. For valence distributions, there is no mixing, 
and their evolution is governed by the gg-kernel alone. The kernels B ab (u,v) 
have the following symmetry 

B ab {u,v) = B ab (v,u) . (66) 

For the parton helicity averaged case, the kernels B ab (v, u) were originally 
obtained in Refs. 37, 38. We present them in the form given in Ref. 16, 

= ^C F (l + 8{u) [v/v]+ + 8(v)[u/u} + -^S(u)S(v)\ , (67) 
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B GQ (u,v) = 


Oi s 

— Of 


( 2+ 


B QG (u,v) = 


-N f (l + 


B 9a {u,v) = 




^4(1 


+ 




'v 2 

V 



S(u)S(v] 
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5{v) 



dv 

v 



+ {u <-> w} 



(69) 



(70) 



Here, /3o = 11— § Nf is the lowest coefficient of the QCD /3-function. Evolution 
kernels for the parton helicity-sensitive case are given in Refs. 39, 40, 

AB qq (u,v) =B qq (u,v) 
AB GQ {u, v) = ^-C F U{u)6{v) - 2\ , 

AB QG (u, v) = — N f (1 - u - v) , 



AB qq (u,v) = B"(u,v) - 12 — N c uv. 

TT 



(71) 
(72) 

(73) 
(74) 



Inserting the operator evolution equation ( pq ) between particular hadronic 
states and parametrizing the matrix elements by appropriate distributions, one 
can get the relevant evolution kernels. In particular, parametrizing nonforward 
matrix element in terms of DDs, one expresses R ab (x,y; £,77; g) in terms of 
B ab (u,v) kernels, e.g., for the gg-kernel one has 



R qq (x, y; f , rj) = ^B qq (y - V x/£, y - fjx/$ . 
In a similar way, one can get the expression for the DGLAP kernel 



P qq {z) 



and for the Brodsky-Lepage kernel 



B qq (l-v- z,v)dv 



(75) 



(76) 



0(y<v) rv 



B qq (v - (y - v)/ri))dv + {y -> y, rj -> 77} . (77) 



6.3 Evolution kernels for SPDs 

The nonforward matrix elements can be also parametrized in terms of SPDs. 
In the case of nonforward parton distributions, the evolution equations have 
the form 



(78) 
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Again, the new kernels Wg (X, Z) can be expressed in terms of the S-kernels, 
e-g-, 



Jo Jo 



v) 5{X -uZ + v{Z - 0) 0(u + v < 1) du dv . (79) 



As we discussed earlier, NFPDs Tq(X) have two components corresponding 
to regions X > £ and A < £. For this reason, one can imagine four different 
possibilities for the kernels W? q (X, Z): 

• both X and Z are larger than £; 

• both A and Z are smaller than £; 

• the original fraction Z is larger than £, but the evolved fraction X is 
smaller than £; 

• Z < ( but A > C 

The last possibility, in fact, is excluded by the delta function in Eq. (|79|). Since 
A = (1 — u — v)Z + v(, we always have A < £ when Z < £. In other words, if 
the initial fraction Z is smaller than £, the evolved fraction A is also smaller 
than £: the parton is trapped in the Z < £ region. 

DGLAP region: Z > (, X > (. Recall, that when A > (, the initial 
parton momentum Xp + is larger than the momentum transfer r = C,p, and 
we can treat the nonforward distribution function Tq{X) as a generalization 
of the usual distribution function /(A) for a somewhat skewed kinematics. 
Hence, we can expect that evolution in the region £<A<l,C<Z<lis 
similar to that generated by the DGLAP equation. In particular, it has the 
basic property that the evolved fraction A is always smaller than the original 
fraction Z . The relevant qq kernel is given by 

l-X/Z 

W q c q (X,Z)\ c < x < z < 1 = ^ J^' Z B m (\l-X/Z-v{l-t/Z)],v\dv. (80) 

Changing the integration variable to w = v(l — (/Z)/(l — X/Z) = v/(l — 
X' /Z'), we obtain the expression in which the arguments of the B-kernels are 
treated in a more symmetric way 

W? q {X, Z)| c <x<^<! = ^— * [ B ab (w (1 - X/Z) ,w(l- X'/Z')) dw , 

zz Jo 

(81) 

where A' = A — £ and Z' = Z — C, are the "returning" partners of the fractions 
X,Z. Moreover, since Z - X = Z' - X', the kernel W? b (X,Z) is given by 
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a function symmetric with respect to the interchange of X, Z with X\Z' . 
This is also true for the gg,qg and gg kernels. Introducing the notation 
Z)\c<x<z<i we present below the explicit expressions 

for the P-kernelsE 



•f" 



5{X - Z) [ dz X 

Jo 1 — z ) 




P? Q {X,Z) = —C 



(82) 
(83) 
(84) 



Pf(A,Z) = — N c { 2 



z 7 



" A) 


r 1 du 


r 1 dv~ 


_2N C 


Jo u 


Jo v _ 



+ S(X-Z) £L- -- - 



(85) 



The formally divergent integrals over u and v provide here the usual "plus"- 
type regularization of the 1/(Z — X) singularities. The prescription following 
from Eq. (|l]) is that combining the 1/(Z — X) and 8{Z — X) terms into 
[Tc_{Z) - Tc_{X)]/(Z - X) in the convolution of P((X,Z) with T C (Z) one 
should change u 1 - X/Z and v -> 1 - X'/Z'. 

In the ( — > limit, the P ( ? b (X, Z) kernels are directly related to the 
DGLAP kernels, 



1 



c (X, Z) - -P QQ (X/Z) , P f ^(X, Z) - —P Qg (X/Z) 
Pf (X,Z) -> §P sQ (X/Z) , Pf (X,Z) - 



z 2 



(86) 



Here one should take into account that the nonforward gluon distribution func- 
tion T^(X) reduces in the limit ( = to Xf g (X) rather than to f g (X). 

In the region Z > (, the evolution is one-sided: the evolved fraction X is 
smaller than the original Z. Furthermore, since if Z < £ then also X < Z, the 
distributions in the X > £ region are not affected by the distributions in the 
X < C regions. Hence, just like in the DGLAP case, information about the 
initial distribution in the Z > ( region is sufficient for calculating its evolution 
in this region. As we will see below, this situation may be contrasted with the 
evolution of distributions in the Z < £ regions: in that case one should know 
the nonforward distribution functions in the whole domain < Z < 1. 

Qualitatively, the evolution in the X 1 Z > Q region proceeds just like in 
the DGLAP evolution: the distributions shift to smaller and smaller values of 
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X. In the DGLAP case, the distributions approach the 6(x) form condensing 
at a single point x = 0. In the nonforward case, the whole region Z < ( works 
like a "black hole" for the partons: after they end up there, they will never 
come back to the X > £ region. 

ERBL region: Z < £, X < £. When £ = 1, the initial momentum 
coincides with the momentum transfer and Tq{X) reduces to a distribution 
amplitude whose evolution is governed by the BL-type kernels, 

W%i x {X,Z) = V ab (X,Z). (87) 

In fact, the nonforward kernels (X, Z) in the Z < £, X < ( region can be 
directly expressed in terms of the BL-type kernels even in the general £^1,0 
case. As explained earlier, if X is in the region X < then the function 
!Fq(X) can be treated as a distribution amplitude ^>q(Y) with Y = X/(. For 
this reason, when both X and Z are smaller than £, we would expect that 
the kernels W£ b (X, Z) must simply reduce to the rescaled BL-type evolution 
kernels V ah (X/(, Z/Q. Indeed, the relation ( |79| ) can be written as 

W^ g (X,Z) = ^^ J B„(u,v)s(j-uj-v{l-^)j 8(u + v<l)dudv. 

(88) 

Comparing this expression with the representation for the V qq (X, Z) kernels, 
we conclude that, in the region where X/( < 1 and Z/C, < 1, the kernel 
Wl q {X, Z) is given by 

W™{X, Z)| < {x , z} < c = - Vio (A7C, Z/C) . (89) 

Transition from Z > £ to X < The BL-type kernels also govern the 
evolution corresponding to transitions from a fraction Z which is larger than 
C to a fraction X which is smaller than £. Indeed, using the (5-function to 
calculate the integral over u in ([79]), we get 

Wl q {X, Z)\ x < ( < z = i B qq ([1 - X/Z - v(l - Q/Z)] ,v\dv, (90) 

which has the same analytic form (^) as the expression for W qq (X, Z) in the 
region X < Z < (. 

As already noted, the evolution jump through the critical fraction £ is 
irreversible: when the parton momentum degrades in the evolution process to 
values smaller than the momentum transfer £p + = r + , further evolution is like 
that for a distribution amplitude: the momentum can decrease or increase up to 
the r + -value but cannot exceed this value. Inside the Z < ( region, the ERBL 
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evolution transforms the ^ (Y) distribution amplitudes into their asymptotic 
forms like YY, YY(Y - Y) for the quarks and (YY) 2 , (YY) 2 (Y - Y) for the 
gluons; a particular form is dictated by the symmetry properties of the relevant 
operators. 



6.4 Asymptotic solutions of evolution equations 

At the leading logarithm (or one loop) lp:el, the solution for QCD evolution 
equations is known in the operator formEj Choosing specific matrix elements 
one can convert the universal solution into four (at least) evolution patterns: 
for usual parton densities ((p| • • • \p) case), distribution amplitudes ((0| . . . \p) 
case), skewed and double parton distributions ({p—r\...\p} case). In the sim- 
plest case of flavor-nonsinglet (valence) functions, the multiplicatively renor- 
malizable operators were originally found in Ref. 5, 

O^ 3 = {zd + ) n i>X a zCl' 2 {z 5 /zd+)i>, (91) 

where C^ 2 (a) are the Gegenbauer polynomials and we use the symbolic no- 
tation (z D /zd+) introduced in Ref. 5, 5=5 -5,9+ =5 + D=d + d- 
In contrast, the usual operators ip\ a z(z D) n tp m ix under renormalization with 

the lower spin operators (zd+) n ~ k Tp\ a z(z D) k ip- In Ref. 5, it was also noted 
that these operators coincide with the free-field conformal tensors. 

Inside the nonforward matrix element, one can substitute (z D) — > (zk) — 
x{zP) and {zd+) — > (zr) = £(zP). Thus, the multiplicative renormalizability 
of 0^ s operators means that the Gegenbauer moments 

C 5 «lM) = C f_Cl' 2 {z/OH NS {z,^)dz (92) 

of the skewed parton distribution H NS (i) have a simple evolution^ 



-•NS 



In /x /A 
ln/i/A 



■y n /0o 



(93) 



where {3q = 11 — |iVy is the lowest coefficienLof the QCD (3-i unction and 7„'s 
are the nonsinglet anomalous dimensions p3a 



7n = Ct 



1 



1 



n+l 



(n + l)(n + 2) 



2 §7 



(94) 



For n — 0, the Gegenbauer moment coincides with the ordinary one and, 
since 70 = 0, the area under the curve remains constant. Other Gegenbauer 
moments decrease as /x increases. 
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Switching from SPDs to DDs, writing the SPD variable x in terms of DD 
variables x = x + a(, and using 

Cl'\x/H + a)=± r ( ( 3 /2) ^ 3/ ^ (2x/0 n - 1 C^ +n -\a) , (95) 

one can express the Gegenbauer moments C n (£,/x) in terms of the combined 
[rr-ordinary <g) a-Gegenbauer] moments of the relevant DDs: 



f J-i J-i+\*\ r(3/2)(n-2fc)! 



x x n - 2k Cl' k 2+n - 2k {a) f NS (x,a\v)da. (96) 

Hence, each x m Cf^ 2+m (a) moment of f NS (x, a; fx) is multiplicatively r- 
malizable and its evolution is governed by the anomalous dimension 7; +r 
In Eq. (p6|), we took into account that the DDs f(x,a) are always even in 
a, which gives an expansion of the Gegenbauer moments in powers of £ 2 . In 
the nonsinglet case, the Gegenbauer moments C„(£,/i) are nonzero for even n 
only. A similar representation can be written for the Gegenbauer moments of 
the singlet quark distributions. In the latter case, the DD f s (x,a) is odd in 
x, and only odd Gegenbauer moments C^(£,/i) do not vanish. 

Another simple case is the evolution of the gluon distributions in pure 
gluodynamics. Then the multiplicatively renormalizable operators with the 
same Lorentz spin n + 1 as in Eq. ( |9~ij ) are 

O a n = z^z^zd+r^G^C^iz D /zd+) Gav ■ (97) 
Due to the symmetry properties of gluon DDs, only Gegenbauer moments 

= C- 1 C C 5 J_\(z/OH G (z,^)dz (98) 



with odd n do not vanish. The Gegenbauer moment can also be written in 
terms of DD: 



[("-l)/2] „l „1-|.t| 

c„ g (^) = e k dx 2«- 2fe - 1 



r(n- 2fc + 3/2) 



^ 7-i r(5/2)(n-2fc-l)! 

x x n - 2k Cl' k 2+n ~ 2k (a) f 9 ( Xl a\^)da. (99) 

Note, that two shifts, n — ► n— 1 and 3/2 — > 5/2, compensate each other. Again, 
each combined x m Cf +m (a) moment of f G {x, a) is multiplicatively renDpnal- 
izable and its evolution is governed by the anomalous dimension 7^^LL3lII 
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Since the Gegenbauer polynomials Cf' 2+m (a) are orthogonal with the 
weighL(l — a 2 ) m+1 , evolution of the x m -moments of DDs in both cases is 
givencil by the formula 

f m {a\^) = J x rn f(x,a\n)dx 

oo 

= (1 - a 2 ) m+1 A ml C? +3/2 (a) [log( M /A)p™+' //3 ° , (100) 

fc=0 

where the coefficients A m i are proportional to x m Cf^ 2+m (a) moments of DDs. 
A similar representation holds in the singlet case, with [log(/i/A)] lm + l ^° 
substituted by a linear combination of terms involving [log(/^/A)] _7m + i ^ and 

[log(///A)] _7m+1 ^° with singlet anomalous dimensions ff^ + i obtained by diag- 
onalizing the coupled quark-gluon evolution equationscil 

Let us consider first two simplified situations. In the quark nonsinglet 
case, the evolution is governed by j^+k a l° ne 

oo 

f» a (a\ /*) = (!- « 2 )" +1 E A nk C n k +3/2 (a) [log( M /A)] 27 ^° . (101) 

fc=0 

Since — while all the anomalous dimensions 7^ with N > 1 are 

negative, only f^ s (a | /z) survives in the asymptotic limit fi — > 00 while all the 
moments f^ s (a\fi) with n > 1 evolve to zero values. Hence, in the formal 
fj, — ► 00 limit, we have 

f NS (x, a I M ^ 00) ~6(x)(l -a 2 ), 

i.e., in each of its variables the limiting function f N (x,a |/x — ► 00) acquires 
the characteristic asymptotic form dictated -by the nature of the variable: 5(x) 
is specific for the distribution functionslpea while the (1 — a 2 )-form_ia. the 
asymptotic shape for the lowest-twist two-body distribution amplitudestm For 
the off-forward distribution of a valence quark q this gives 

H wd n(.x,Z\p^ o) = 1(1 -&/?)/£. 

Another example is the evolution of the gluon distribution in pure gluo- 
dynamics which is governed by 7^_ fe with /3q = UN c /3. Note that the lowest 
local operator in this case corresponds to n — 1. Furthermore, in pure gluody- 
namics, jf 9 vanishes while "fff < if N > 1. This means that in the /1 — > 00 
limit we have 

xf 9 (x,a\fx^ 00) = ^5(x)(l- a 2 ) 2 
lb 
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for the double distribution which results in 

for the off-forward distribution. In the formulas above, the total momentum 
carried by the gluons (in pure gluodynamics!) was normalized to unity. 

In QCD, we should take into account the effects due to quark gluon mix- 
ing. The diagonalization procedure gives two multiplicatively renormalizable 
combinations 

P nk = fnk + c nkfnk 

where (omitting the nk indices) 



= 7 3g - y?q ± ^/frflfl - 7 -?9)2 + ^gg^g 

Their evolution is governed by the anomalous dimensions 

,± 



\k = \ (7" + l qq ± V(l" - l qq ) 2 + 47^7^ . (104) 



In particular, 7^ = and af = 1 which means that the sum ff + ff does 
not evolve: the total momentum carried by the partons is conserved. Another 
multiplicatively renormalizable combination involving ff and ff is 



Jw 4N f Jw ' 



It vanishes in the \i — > 00 limit, and we have 

- »> ~ IcTTW, ' '><•<■» - ~> - T^TW, ■ < 105 > 

Since all the combinations P nk with n + k > 2 vanish in the /x — > 00 limit, we 
have 

xf{x, a I M - 00) - H 4 ^ F ^ - a 2 ) 2 , (106) 

a I M - 00) - ^ JcTTWf 5{X){1 ~ a2f ' (1 ° 7) 

or 

f\x, a I M - 00) 4C j^ <y'(x)(l - a 2 ) 2 • (108) 

In terms of off-forward distributions this is equivalent to 

f 1 ^ 00) - H (i _ , V a 2 /e , (109) 
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H*(x,t | /i - oo) - ^ *(! - ^/e 2 ) 2 /C 2 • (HO) 

Note, that in the /i — > oo limit both the functions H q ' 9 (x, £ \ fi) and their 
derivatives (d/dx)H q ' 9 (x,£ : | /x) vanish at the border-point x = £. 



6. 5 Reconstructing SPDs from usual parton densities 

The anomalous dimensions 7„ increase with raising n, and, hence, the mth x- 
moment of fix, a; fi) is asymptotically dominated by the a-profile (1 — a 2 ) m+1 . 
Such a correlation between x- and a-dependences of f(x, a; /i) is not something 
exotic. Take a DD which is constant in its support region. Then its x m - 
moment behaves like (1 — |a|) m+1 , i.e., the width of the a profile decreases 
with increasing n. This result is easy to understand: due to the spectral 
condition \a\ < 1 — \x\, the x m moments with larger m are dominated by 
regions which are narrower in the a-direction. 

These observations suggests to try a model in which the moments / to (q; /x) 
have the asymptotic (1 — a 2 ) m+1 profile even at non-asymptotic \i. This is 
equivalent to assuming that all the combined moments x rn Cf^ 2+m (a) with 
I > vanish. Note that this assumption is stable with respect to PQCD 
evolution. Since integrating f m (a; fi) over a one should get the moments 
f m (fi) of the forward density f(x; /z), the DD moments / m (oi; /j) in this model 
are given by 

fm{a;jj) = Pm+l(<x)fm(jt) ( m ) 

where p m +i(a) is the normalized profile function (cf. Eq.(p4|)). In the explicit 
form: 

r(m + 5/2) n +1 << 



x m f(x, a;»)dx= r(1 ) 2) m + 1 )! (1 " * } J_ x ' (H2) 

In this relation, all the dependence on a can be trivially shifted to the left-hand 
side of this equation, and we immediately see that fix, a ; /i) in this model is 
a function of x/(l — a 2 ), 

fix, a; fi) = F(x/(1 - a 2 ); (i) 6>(0 < x/(l - a 2 ) < 1) . (113) 

A direct relation between /(z, /i) and F{u; /x) can be easily obtained using the 
basic fact that integrating fix, a; fi) over a one should get the forward density 
f(z,fi); e.g., for positive z we have 

m=z[ 1 ^ { ^— du. aw) 
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This relation has the structure of the Abel equation. Solving it for F(u) we 
get 

u*' 2 f 1 [f(z)/z]' 



F(u) = / LJ ; " 1 dz. (115) 

n Ju \/z-u 

Thus, in this model, knowing the forward density f(z) one can calculate the 
double distribution function f(x,a) = F(x/(1 — a 2 )). 

Note, however, that the model derived above violates the DD support 
condition |x| + |a| < 1: the restriction |x| < 1 — a 2 defines a larger area. Hence, 
the model is only applicable in a situation when the difference between two 
spectral conditions can be neglected. A practically important case is the shape 
of H{x,£) for small £. Indeed, calculating H(x,£) for small £ one integrates 
the relevant DDs f(x) over practically vertical lines. If x is also small, both 
the correct \a\ < 1 — \x\ and model ct 2 < 1 — |x| conditions can be substituted 
by \at\ < 1. Now, if x 3> £, a slight deviation of the integration line from the 
vertical direction can be neglected and H(x,^) can be approximated by the 
forward limit f(x). 

Specifying the ansatz for /(z), one can get an explicit expression for the 
model DD by calculating F(u) from Eq. ( |115| ). However, in the simplest case 
when f(x) — Ax~ a for small x, the result is evident without any calculation: 
the DD f(x, a) which is a function of the ratio x/(l — a 2 ) and reduces to Ax~ a 
after an integration over a must be given by 

f(x,a) = p a (a)f(x) , 

where p a {ct) is the normalized profile function of Eq. (p4|), 

f(x,a)=A F i"t 5/2) A l~a 2 ) a x- a . (116) 
Jy ' ' r(l/2)T(a + 2) v ; v ' 

This DD is a particular case of the general factorized ansatz f(x, a) = 
Pn{ot)f(x) considered in the previous section. Its most nontrivial feature is 
the correlation n = a between the profile function parameter n and the power 
a characterizing the small- a: behavior of the forward distribution. 

Knowing the DDs, the relevant SPDs H a (x,ti) can be obtained in the 
standard way from f a {x, a) for quarks and from xf q {x, a) in the case of gluons. 
In particular, the SPD enhancement factor 1Z(Q for small £ in this model is 
given by 

*C(0 + (n7) 



for quarks and by 



fi(0 r(a + 2)r(a + l) 
J?(C) r(2a + 2) 

~ r(a + 3)r(a+i) 



(118) 
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for gluons. 

The use of the asymptotic profiles for DD moments is the basic 

assumption of the model described above. However, if one is interested in 
SPDs for small £, the impact of deviations of f n (ot) from the asymptotic profile 
is suppressed. Even if the higher harmonics are present in /„(«), i.e., if the 
x n-2k(j3/2+n 2k ^ moments of f(x,a) are nonzero for k > 1 values, their 
contribution into the Gegenbauer moments C„ is strongly suppressed by 
£ 2fc factors (see Eq. (|96|)). Hence, for small £, the shape of H(x,£) for a wide 
variety of model a-profilcs is very close to that based on the asymptotic profile 
model. 

Absence of higher harmonics in fn{ct) is equivalent to absence of the £- 
dependence in the Gegenbauer moments C n (£, fi). The assumption that the 
moments C n (£, /i) do not depend on £ was the starting point for the model of 
SPDs H(x,£) constructed in Ref. 43. Though the formalism of DDs was not 
used there, both approaches lead to identical results: the final result of Ref. 43 
has the form of a DD representation for H(x,£,). 



7 DVCS amplitude at leading twist and beyond 

7.1 Twist-2 DVCS amplitude for the nucleon 

Using the parametrization for the matrix elements of the quark operator, we 
can easily write the leading twist contributionOO to the DVCS amplitude 



dx 



1 x - £ + iO 
1 



1 



-5^ + ^7y(W + ^V M ) 



x\u(p')q'u{p)H a (x^;t) + —u(p')(q'f - f^)u{p)E a {x,i;t) 
(q'r) 



+ie ^_^_ j u(p')q' l5 u(p)H a (x,Z;t) 



2M 



u(p')-y 5 u(p)E a (x,£,;t) 



(119) 



Note, that the functions H a (x, £; f), E a {x, £; t) parametrizing the matrix ele- 
ment of the O a operator are odd in x while the distributions H a {x,^\t) and 
E a {x^;t) related to 0^ a term are even in x. Alternatively, one can use the 
combinations 

^[i/(i-e + *o)±i/(s + e+to)] 

in which the contributions of the s- and it-channel handbag diagrams are ex- 
plicitly separated. 
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Thus, the skewed parton distributions appear in the DVCS amplitude in 
an integrated form. Note that the relevant integrals 

r 1 dr 
h a (Z',t)= / H a (x,Z;t) * (120) 
J— l x — ^ -\- l\J 

have both real and imaginary parts. The latter are given by the values of the 
relevant SPDs at the border point x = £ 

hnh a {t;t) = -*H a (t,&t) . (121) 

For a fixed t, the "effective" SPD H a (£, £;t) is a function of the Bjorken ratio 
xsj = 2£/(l + £), just like DIS structure functions. A linear combination of 
the effective (or border-point) SPDs -ff a (£, £; tLis directly accessible through 
the measurement of the single-spin asymmetryuil Another function of x Bj cor- 
responds to the real part of h a (^;t). It is given by the principle value integral 

r 1 dx 

Refc„(6t) = P/ H a (x,Z;t)^—. (122) 

The real part of the DVCS amplitude can be accessed through the measurement 
of the lepton charge asymmetrytil 



In Eq. (119), the final photon momentum q' is used as a natural light- 
cone 4-vector specifying the "minus" direction. In this form, the amplitude 
T^ v (P,r,q') exactly satisfies the transversality condition T fJ,t, (P,r,q')q' u = 
with respect to the final photon momentum. However, the convolution 
T fl, ' J (P,r,q')q ll is proportional to r p — P^{rq') / {Pq') = Aj_, the transverse 
component of the momentum transfer A = r±. Hence, the accuracy of the 
twist-2 approximation is not sufficient to satisfy the transversality condition 
T MI/ (P, r, q')q^, = 0. Guichon and Vanderhaeghen E3 (GV) proposed to add a 
non-leading O(A) term producing the expression 

= T?:_ 2 - ^( ?A T&_ 2 ) (123) 

which satisfies both q^T^ = and q' u T^ = 0. 

It is important to note that the use of the GV prescription changes the 
{/j v\ symmetry structure of the DVCS amplitude. In particular, the GV 
expression constructed from the {fi «-> v\ symmetric part of T^ w _ 2 satisfies the 
transversality conditions but it is not symmetric in \xv anymore. It is easy to 
see that this is a common feature. Indeed, the transversality conditions written 
in symmetric variables Q = (q + q')/2, r = q' — q and P = (jp + p')/2 convert 
into two relations 
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connecting the symmetric T^ v y = (!),„ + T vll )/2 and antisymmetric 2~L„i = 
(T^ — T v ^)j2 parts of T MI/ . In the r = forward limit, the two relations 
decouple to give the DIS transversality conditions q^Tr^x = 0, gTu,] = 0. 

The GV prescription was supported by several groups 46 " 48 ' 24 who derived 
this term in a regular way as a kinematical twist-3 contribution. Note, that 
the twist-3 quark-gluon operators qGq are dynamically independent from the 
twist— 2 qq ones. Hence, to get a gauge invariant extension of the twist-2 
contribution, it is sufficient to retain only the part of the twist-3 SPD's induced 
by the twist-2 distributions, i.e., the Wandzura-Wilczek (WW) type terms. 

A very convenient way to analyze the DVCS amplitude bevand the leading- 
twist level is provided by the approach of Balitsky and Braunt3 (see, however, 
Refs. 48-51 where other versions of the light cone analysis are used). We 
combine it with the formalism of double distributionsEJH which provides a 
simple way of deriving relations between SPD's describing the kinematical 
twist-3 effects and the basic twist-2 distributions. 



7.2 Twist decomposition 

The nonlocal operators O a , 0^ a in Eq. (g) do not have a definite twist. The 
twist-2 part of these operators is defined by formally Taylor-expanding the 
nonlocal operators in the relative coordinate z and retaining only the totally 
symmetric traceless parts of the coefficients in the expansion: 



| Z M1 



oo 

$(X z/2hMX + z/2)] twist ~ 2 = £ 1 

71=0 

tp(X) 7{CT B M1 • ■ • D„ n} -traces i/)(X) , (125) 

and similarly for the operator with 7^75 (cf., e.g., Ref. 41). The sjanmetrization 
can be carried out directly at the level of non-local operatorsEj Indeed, the 
part of the nonlocal operator corresponding to totally symmetric local tensor 
operators is projected out by 

[${X - z/2) la tP(X + z/2)] Bym = [ dt$(X -tz/2)zi>(X + tz/2) . 

uz a J Q 

(126) 

The subtraction of traces in the local operators implies that the twist-2 string 
operator contracted with z a should satisfy the d'Alembert equation with re- 
spect to z, 

a z [i>{X - tz/2)z^{X + tz/2)] twist ~ 2 = 0. (127) 



In the center-of-mass X and relative z coordinates, the transversality con- 
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ditions (124) are 



An w (z|x) = i^|-n M (z|x) , (128) 
jL-n W] { z \x) = i JLn {|i „ } (*|x) . (129) 

Consider the part of the current product given by Eq. (p|) with the nonlocal 

rom Eq 



operators replaced by their twist- 2 parts. From Eq. (127) and 

d z 



dz p 2tt 2 z 4 



it follows that 



Fl 8 

Tjtwist — 2 ri yrtwist — 2 ri / -i Qfl\ 

Since forward matrix elements are zero for all total derivative operators, this 
guarantees the transversality of the twist-2 contribution in the case of deep 
inelastic scattering. In the non-forward case, we have 

(d/axjiifift- 2 * o, (d/dx^™*- 2 ± , 

and ( |129| ) is violated. The non- transverse terms in the twist-2 contribution 
can only be compensated by contributions from operators of higher twist. In 
fact, the necessary operators are contained in the part of the string operator 
which was dropped in taking thft-twist-2 part. Incorporating QCD equations 
of motion, it is possible to showed that the twist > 2 part involves the total 
derivatives of nonlocal operators 

^(-z/2M(z/2) - [^(-z/2) 7t ^(z/2)] sym 

= 2 eaipKH dX~ J Q dtt $(- tz / 2 )W^{tz/2) + ... . (131) 

The ellipses stand for quark-gluon operators (we do not write them explicitly 
since they are not needed to restore transversality of the twist-2 contribution, 
but, in principle they can be kept). The relation for the operator containing 
7 a 75 is obtained by changing j a — > 7^75, 7k75 — * 7«- 

Note, that the operators appearing under the total derivative on the right 



hand side of Eq. ( |131| ) and its 7q75 analog are still the full string operators with 
no definite twist. Hence, one can decompose them into a symmetric (i.e., twist- 
2) part and total derivatives, and so on; thus expressing the original string 
operator as the sum of its symmetric part and an infinite series of arbitrary 



44 



Handbook of QCD / Volume 2 



order total derivatives o£ symmetric operators. This series can be summed 
up in a closed formE3ejE-3 Up to operators whose matrix elements give 0(t) 



ivcs_oi_s 

contributions to the Compton amplitude, the result isEa 



-2/2)7^/2) = / 




Jo 


iv 

+ -sm 


iv f d 
.2 \ Z dX 


i 


d d 





dv < cos 



iv 
~2 



dX 



d 

dz a 



d 



l 

dv 



ip(-vz/2)zip(vz/2) 



du cos 



iv. 

T 



dx 



ip(-vz/2)z'y 5 'ip{vz/2) + 



(132) 



(see also Refs. 48, 49). An analogous formula applies to the operators with 
7<t —> 7<t75; one should just replace z 275, 275 — * z. 



7.3 Parametrization of nonforward matrix elements 

Double distributions. To get the amplitude for deeply virtual Compton scat- 
tering off a hadronic target we need parametrizations of the hadronic matrix 
elements of the uncontracted twist-2 strin g op erators O a ,0^ G appearing in 
Eq. (||). We will derive them from Eq. (132). For simplicity, we consider 
here one quark flavor and the pion target, which has zero spin and practically 
vanishing mass. In this case, the matrix element of the contracted axial oper- 
ator z a 0^ a (z I 0) (parametrized in the forward limit by the polarized parton 
density) is identically zero. Thus we need only the parametrization for the 
matrix element (P — r/2 | 0(z \ 0) | P + r/2) of the contracted vector operator 
0(z I 0) = z"O a {z I 0). With respect to z, it can be regarded as a function of 
three invariants (Pz), (rz) and z 2 . For dimensional reasons, the dependence on 
z 2 is through the combinations tz 2 and P 2 z 2 only. We are going to drop 0(t) 
and 0(P 2 ) terms in the Compton amplitude, so we may ignore the dependence 
on z 2 and treat this matrix element as a function of just two variables (Pz) and 
(rz). Incorporating the spectral properties of nonforward matrix elementspa 
we write the plane wave expansion in the form 



(P -r/2\O(z\0)\P + r/2) 



2(Pz) 



+ (rz) 



,l-|x| 

dx I e~ l{kz) f(x,a)da 
-j j— 



-<a(r*)/2 £,( a ) da ) 



(133) 



where k = xP + ar/2, f(x,a) is the double distribution (DD) and D(a) 
is the Polyakov- Weiss (PW) distribution amplitudeEil absorbing the (Pz)- 
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independent terms. From this parametrization, we can obtain matrix elements 
of original uncontracted operators, ( |132| ), including the kinematical twist-3 
contributions. We consider first the part coming from the double distribu- 
tion term in Eq. ( 133 ); the contributions from the PW-term will be included 
separately. In matrix elements, the total derivative turns into the momentum 
transfer, id/dX a — > r a — 2£P cr +A (T . Similarly, we write k = (x+£_a)P+aA/2. 
This gives 



1 rl-\x 



\{P-r/2\0 !T {z\Q))\P + r/2)= I dx I daf(x,a) 

1 J-l J-l+\x\ 

x jp CT e~ l(£+ « Q)(P2) - M(Az)/2 + ~ [A a {Pz) - P a {Az)] 
x J dvve -iv(x+t a )(Pz)-iv a (Az)/2 [ sin ( 5 ( rz )/ 2 ) _ iacos(v(rz)/2)] J . (134) 

Skewed distributions. Expanding exp[— ia(Az)/2] — 1 — ia(Az)/2 + . . . and 
keeping only terms up to those linear in the transverse momentum A we get 
an expression^] in which the spectral parameter x appears in the exponential 
factors only in the combination x = x+£a. Thus, we can introduce two skewed 
parton distributions, 

i ,i-|x| r i 

dx / da 8(x — x — £a) f(x, a) < . (135) 
-l J-i+\x\ [ a 

Note, that in our case the DD f(x,a) is even in a and odd in x. As a result, 
the functions H and A satisfy the symmetry relations 

H(x,S)=-H(-x,0 , H(x,0=H(x,-0 , 
A(x,0 = A(-x,0 , A(x,0 = -A(x,-0 ■ (136) 

Furthermore, because of the antisymmetry of the combination af(x,a) with 
respect both to x and a we have 




/ dxA{x,£) = 0. 
Jo 



(137) 



Hence, the distribution A(x,£) cannot be a positive-definite function on < 
x < 1. 

Uniting the cosine and sine functions with the overall exponential factor 
e -ivx(Pz) Qne g e j. g vx j_ combinations. Using (136), one can arrange that 

''Because of this truncation, the A^A,, terms in the expression f or th e amplitude T^ v 
will be lost. If needed, they can be kept; see the discussion after Eq. (1143). 
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only vx + v£ would appear, 

I (P - r/2 | | 0) | P + r/2) = P CT / dx e -" (pz) 



-- [A a (Pz) - P CT (Az)] y dx [if (i, - A(x, 0] 
x / du v cos[(?;a: + u£)(Pz)] . 



(138) 



In a similar fashion, we get parametrization for the matrix element of the axial 
string operator ( |l32] ), 

1 i f 1 

- (P - r/2 | 5a (z | 0) | P + r/2) = - e aa/3l z a A 0Pl J dx [H{x, - A(x, £)] 



x / dv-u sin[(wx + v£)(Pz)] . 



(139) 



Note that it is expressed in terms of the same skewed distributions H(x,£) 
and A(x,£) which, in turn, are determined by the original double distribution 
f(x, a), see Eq. (135). 



7.4 DVCS amplitude for pion target 



DD-generated contribution. Substituting the parametrizations ( J138| ) and (139) 
into Eq. (|^) and performing the Fourier integral over the separation z one 
obtains the Compton amplitude 



T — — 



P^Qu + Q^Pu - g^(PQ) + 2£P M P„ + ^-P v -P^ 



dx 



(PQ) 

x - £ + tO (PQ) 



1 



dx R(x, £) 



1 dy (Q„ + 3^)A v 



£ + vx + v£ — iO 



(PQ) J-i 



dxR(x,^) / dv 



A^(Q V + (P V ) 



(140) 



-£ + vx + vt; + i0 ' 

where R(x, £) is a new SPD describing the kinematical twist-3 contributions, 

dH(x,0 dA(x,0 



R(x,0 



dx 



dx 



(141) 



All three terms in Eq. (140) are individually transverse up to terms of order 
i,P 2 . 

Singularities. The first term is the twist-2 part with the tensor structure 
corrected exactly as suggested by Guichon and Vanderhaeghenc3 The integral 
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over x exists if H(x,£) is continuous at x — £, which is the case for SPD's 
derived from the DD's that are less singular than 1/x 2 for x = and are con- 
tinuous otherwise (see Ref. 23). In particular, continuous SPD's were obtained 
in model calculations of SPD's at a low scale in the instanton vacuumEa The 
second term contributes only to the helicity amplitude for a longitudinally 
polarized initial photon. The parameter integral over v gives the function 
[ln(x + £ — iO) — ln(2£ — iO))/ (x — £) which is regular at x = £ and has a loga- 
rithmic singularity at x — — £. The integral over x exists if R(x, £) is bounded 
at x = — £, which again is the case in the DD-based models described in Ref. 



23. The third term of Eq. ( 140 ) corresponds to the transverse polarization of 
the initial photon. In this case, one faces the integrand l/[v(x — £) + iO] which 
produces dv/v divergence for the u-integral at the lower limit. One may hope 
to get a finite result only if the integral 

1(0 ee f dx R{x f (142) 
J -i x-E, + iQ 

vanis hes. From the definition of the skewed distributions H(x,£) and A(x, £) 
( |135| ) it follows that 

dA(x,Q _ dH(x,Q 
dx d£ 

Hence, one can substitute R(x, £) by the combination dH{x,£)/dx + 
dH{xUW£ (see Refs. 48, 49, 24). Integrating the dH{x,£)/dx term by 
parts gives 

x — 4 + iQ 



i.e., the £ derivative of the twist-2 contribution. In general, the latter has a 
nontrivial ^-dependent form determined by the shape of SPDs (see, however, 
the discussion of the PW contribution below). Hence, the twist-3 part of the 
tensor amplitude T^ v diverges in case of the transverse polarization of the ini- 
tial photonoEJ However, it is easy to see that the relevant tensor structure 
&n{Qv +£Pv) is just a truncated version of the exactly gauge invariant combi- 
nation A^q' v which has zero projection onto the polarization vector e' v of the 
final real photon: (e'q 1 ) = 0. 

The structur e A u q' v is obtained if one uses the original full form of the DD 



parametrization (134). It appears from the term with the exponential factor of 
the argument — i[v(x+£a)(Pz)+v£(Pz) + (va+v)(Az)/2] which is obtained by 
combining the sine/cosine functions and the exponential factor in the second 
term of Eq. ( |134[ ). In the Compton amplitude, it gives rise to a contribution 
in which the argument of the quark propagator is q + v(x + £a)P + vt;P+ (va + 
v)A/2. Since (AQ), (AP) and A 2 are negligible, the denominator factors in 
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Eq. ( |140p remain unchanged. In numerators, representing (to + v)A/2 as 
[1 — (1 — a)v]A/2, we observe that A^(Q L , + ^P v ) converts into the A fl (Q L/ + 
£P, y + A„/2) — A^q' v term plus a vA^A v type contribution corresponding 
to a new SPD built from the (1 - a) 2 f(x,a) DD (cf. fll~35|)). Due to the 
extra v factor, the v-integral for the latter contribution is finite. Hence, for the 
physical DVCS amplitude, we find no evidence for factorization breaking in 
the kinematical twist-3 contributions, both in their 1/ — q 2 and l/q 2 terms. 
It is quite possible that factorization breaks down at the l/q 2 level, but one 
needs to analyze 0(z 2 ) suppressed terms (i.e., twist-4 contributions) to see if 
it really happens. 

\L±upe representation. These results can be expressed in another 
3EJ by introducing new skewed distributions related to R(x,£/) via 
the integral transformation similar to that used by Wandzura and WilczekE^I 
Treating the combination xv + in (140) as a new variable we define 



for 



R(y,O d V / S(yv + v£-x)dv 



(x>0 



1 gy 



dy - 0(x < 



-i y-i 



dy. (144) 



In terms of this transform, the matrix element of the vector operator (138) can 
be expressed as 



1 f 1 

-(P-r/2\O (7 (z\0)\P + r/2) = dx 

(Az) 



- l -P a {Az)A{x,i) + \ 



A a -P a 



(Pz) 



Rw(x,£) - Rw{-x,£) 



.(145) 



Note that only the odd part of Rw{x,£) contributes here. In case of the axial 
operator (139) 



-(P-r/2\O 5a (z\0)\P + r/2) 



4 £ffQ/37 (Pz) 



A P- t 



dxe- lx(Pz) 



Rw{x, + Rw(-x, £) 



(146) 



only the even part of Ryy(x,£) appears. The part of the Compton amplitude 
(140) containing R(x, £) can be written in terms of this new function as 



(QP) 



A»(Q V + £P V ) (Qf + 3^P")A V 



£-i0 



Rw(x, £) dx 



(147) 



The integrals with l/(x ± £ =F iO) converge only if the function Ry/(x,£) 
is continuous for x — ±£. According to Eq. (U44), £) is given by the 
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integral of R(y,£)/(y — £) from x to 1 if x > £ and from a; to —1 if .t < £. 
Evidently, a; = — £ is not a special point in the integral transformation (144), 
hence the function R\y(x, £) is continuous at x = — £. However, it is extremely 
unlikely that the limiting values approached by Rw{ x > for x = £ from below 
and from above do coincide. Indeed, JJm, difference of the two limits can be 
written as the principal value integral 



Rw(t + 0,0- Rwti - 0, = P \ dy , 



(148) 



which can be converted into the ^-derivative of the real part of the twist-2 
contribution. This means that the singularity, which was obtained as a straight 
divergence of the dv/v integral, in the WW- type approach appears due to an 
unavoidable discontinuity of the Rw(x,£) transform at x = £. 

Contribution from the PW term. The contribution of the PW term to the 
vector operator 



^(P-r/2\O a (0\z)\P + r/2) PW term 



= [ e- ia{rz)/2 D{a)da (149) 
2 J-i 



has a simple structure corresponding to a parton picture in which the partons 
carry the fractions (1 ± ot)/2 of the momentum transfer r. Since only one mo- 
mentum r is involved, this term can contribute only to the totally symmetric 
part of the vector string operator: it "decouples" in the reduction relations 
(131). In particular, the PW term does not contribute to the second contri- 
bution in Eq. (132) which is generated by decomposition of the axial string 
operator: both derivatives, with respect to X and z, give rise to the momentum 
transfer r, whence the contraction with the e-tensor in (132) gives zero. Thus, 
the PW-contribution should be transverse by itself. Indeed, a straightforward 
calculation gives 



T 



PW 



2 
(rq) 



1 D(a) 



da , (150) 



liu |PW 



= 0. Hence, this term can 



which evidently satisfies g^T^lpw = 0, 
be treated as a separate contribution. 

Alternatively, one may include it into the basic SPD H(x,£,) and all 
SPD's derived from H(x,£). Specifically, for £ > 0, the PW contribution 
to H(x,§ is D{x/i) 6{\x\ < Z)fB it contributes (£ - x)D'{x/(,) 6(\x\ < £)/£ 2 
[where D'(a) = (d/da)D(a)] to R(x, £); furthermore, the PW contri butio n to 
Rw( x , Q is D{x/0) 6{\x\ < £) /£■ Inserting these functions into Eqs. (140) and 
(147) one rederives Eq. ( |150[ ). One can also observe that the PW term gives 
zero contribution into /(£), Eq. (142). 
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8 Real Compton scattering 

8.1 Compton amplitudes and light-cone dominance 

The Compton scattering in its various versions provides a unique tool for study- 
ing hadronic structure. The Compton amplitude probes the hadrons through a 
coupling of two electromagnetic currents and in this aspect it can be considered 
as a generalization of hadronic form factors. In QCD, the photons interact with 
the quarks of a hadron through a vertex which, in the lowest approximation, 
has a pointlike structure. However, in the soft regime, strong interactions pro- 
duce large corrections uncalculable within the perturbative QCD framework. 
To take advantage of the basic pointlike structure of the photon-quark cou- 
pling and the asymptotic freedom feature of QCD, one should choose a specific 
kinematics in which the behavior of the relevant amplitude is dominated by 
short (or, being more precise, lightlike) distances. The general feature of all 
such types of kinematics is the presence of a large momentum transfer. For 
Compton amplitudes, there are several situations when large momentum trans- 
fer induces dominance of configurations involving lightlikc distances: 

i) both photons are far off-shell and have equal spacelike virtuality: virtual for- 
ward Compton amplitude, its imaginary part determines structure functions 
of deep inelastic scattering (DIS); 

ii) initial photon is highly virtual, the final one is real and the momentum 
transfer to the hadron is small: deeply virtual Compton scattering (DVCS) 
amplitude; 

Hi) both photons are real but the momentum transfer is large: wide-angle 
Compton scattering (WACS) amplitude. 

The first two cases were discussed in the previous sections. As argued in 
Ref. 25, at accessible momentum transfers \t\ < 10 GeV 2 , the WACS ampli- 
tude is also dominated by handbag diagrams, just like in DIS and DVCS. In 
the most general case, the nonperturbative part of the handbag contribution 
is described by double distributions (DDs) f(x, a; t), g(x, a; t), etc., which, as 
discussed earlier, can be related to usual parton densities f(x), Af(x) and 
nucleon form factors like Fi(t), G^t). Among the arguments of DDs, x is the 
fraction of the initial hadron momentum carried by the active parton and y is 
the fraction of the momentum transfer r. The description of the WACS am- 
plitude simplifies when one can neglect the a-dependence of the hard part and 
integrate out the a-dependence of the double distributions. In that .case, the 
long-distance dynamics is described by nonforward parton densities cB (NDs) 
T{x\ t), Q{x; t), etc. The latter can be interpreted as the usual parton densities 
fix) supplemented by a form factor type t-dependence. A simple model for the 
relevant NDs was proposed in Ref. 25. It both satisfies the relation between 
T(x\ t) and the usual parton densities f(x) and produces a good description of 
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the F\{t) form factor up to t ~ —10 GeV . This model was used in calculate 
the WACS amplitude in a rather close agreement to existing datal 



>3 



8.2 Modeling nonforward densities 

Let us apply the DD formalism to the large-i real Compton scattering. Since 
the initial photon is also real, q 2 = (and hence xsj — 0), it is natural to 
expect that the nonperturbative functions which appear in WACS correspond 
to the £ = limit of the skewed parton distributions. In the £ = limit, the 
SPDs reduce to the nonforward parton densities 

/l — \x\ 
f a (x,a;t)da. (151) 
-1+1*1 

Note that NDs depend on "only two" variables x and t, with this dependence 
constrained by reduction formulas 

j*( x . t = 0) = f a ( x ) ; VeJ [T a {x-t)-r a {x-t)]dx = F 1 {t). (152) 

Furthermore, it is possible to interpret the nonforward densities in terms 
of the light-cone wave functions. Consider for simplicity a two-body bound 
state whose lowest Fock component is described by a light cone wave function 
^(x, k±). Choosing a frame where the momentum transfer r is purely trans- 
verse r = rj_, we can write the two-body contribution into the form factor 



asi 



F t>,r„ /.-../„ {/i= j ,/. r / **(x,k ± +xr ± )*(x,k ± )^±, (153) 



o 

where x = 1 — x. Comparing this expression with the reduction formula ( |152[ ), 
we conclude that 

■$*(x,k ± +xr ± )V(x,k ± ) i (154) 

167r J 

is the two-body contribution into the nonforward parton density T{x, t). As- 
suming a Gaussian dependence on the transverse momentum k± (cf. Ref. 54) 

■^{x,k ± ) = ^{x)e- k ^' 2xSX \ (155) 

we get 

■pitwo-body) t ) = j (two-body) ^ £ xt/4x\ 2 ^ ^gg) 
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where 

f (two-body)^ = $2 = ^ (i6) = 

lOTT 2 

is the two-body part of the relevant parton density. Within the light-cone 
approach, to get the total result for either usual f(x) or nonforward parton 
densities !F(x,t), one should add the contributions due to higher Fock com- 
ponents. These contributions are not small, e.g., the valence du contribution 
into the normalization of the ir + form factor at t = is less than 25% (see 
Ref. 54). In the absence of a formalism providing explicit expressions f or an 



infinite tower of light-cone wave functions one can choose to treat Eq. (156) 
as a guide for fixing interplay between the t and x dependences of NDs and 
propose to model them by 

F a ( X ,t) = f a (x)e St ^ xX2 = I e -(ki + (^+*r ± f)/2 xS X> d 2 k± ( 15 g) 



nxx 



The functions f a (x) here are the usual parton densities. They can be taken 
from existing parametrizations like GRV, MRS, CTEQ, etc. In the t = 
limit (recall that t is n egati ve) this model, by construction, s atisf ies the first 



of reduction formulas (152). Within the Gaussian ansatz (158), the basic 
scale A specifies the average transverse momentum carried by the quarks. In 
particular, for valence quarks 

(k 2 ± ) a = ^-j\xr a al (x)dx, (159) 
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where N u = 2, — 1 are the numbers of the valence a-quarks in the proton. 

The magnitude of A can be fixed using the second reduction formula in 
( |152| ) relating nonforward densities T a {x,t) to the Fi(t) form factor. The 
following simple expressions for the valence distributions 

fZ a \x) = 1.89a;- - 4 (l-x) 3 - 5 (l + 6a;), (160) 

fi al (x) = 0.54a;- - 6 (l-s) 4 - 2 (l + 8a;). (161) 
closely reproduce the relevant curves given by the GRV parametrizationEl at 



2.5 




-t (GeV") 

Figure 10: Ratio Ff(t)/D(t) of the Ff(t) form factor to the dipolc fit D(t) = 1/(1 - 
i/0.71 GeV 2 ) 2 . Curve is based on Eq. (47) with A 2 = 0.7 GeV 2 . Experimental data are 

taken from Ref. 56. 

a low normalization point Q 2 ~ 1 GeV 2 . The best agreement between the 
model ^ 

FZ° lt (t) = [ [e u fl al {x) + e d fr\x)] e^'^ dx (162) 
Jo 

and experimental dataS in the moderately large t region 1 GeV 2 < \t\ < 10 
GeV 2 is reached for A 2 = 0.7 GeV 2 (see Fig. |o|). This value gives a reasonable 
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magnitude 

{k\) u = (290 MeV) 2 , (kl) d = (250 MeV) 2 (163) 

for the average transverse momentum of the valence u and d quarks in the 
proton. 

Similarly, building a model for the parton helicity sensitive NDs G a (x, t) = 
H(x,£, — 0;i) one can take their t = shape from existing parametrizations 
for spin-dependent parton distributions Af a (x) and then fix the relevant A 
parameter by fitting the G^i) form factor. The case of hadron spin- flip dis- 
tributions JC a {x,t) = E(x,£ = 0;f) and T a (x;t) = E(x,£ = 0;t) is more 
complicated since the distributions e a (x), e a (x) are unknown. 

At t — 0, the model by construction gives a correct normalization 
Ff(t = 0) = 1 _£or the form factor. Moreover, the curve is within 5% from 
the data pointsE3 for 1 GeV 2 < — t < 6 GeV 2 and does not deviate from 
them by more than 10% up to 9 GeV 2 . Modeling the i-dependence by a 
more complicated formula (e.g., assuming a slower decrease at large t, and/or 
choosing different A's for u and d quarks and/or splitting NDs into several 
components with different A's, etc., (see Ref. 30 for an example of such an 
attempt) or changing the shape of parton densities f a {x) one can improve the 
quality of the fit and extend agreement with the data to higher t. However, 
the very fact that a reasonable description of the F% (t) data in a wide region 
1 GeV 2 < \t\ < 10 GeV 2 was obtained by fixing just a single parameter A re- 
flecting the proton size is an evidence that the model correctly catches the 
gross features of the underlying physics. 

Since the model implies a Gaussian dependence on the transverse momen- 
tum, it includes only what is usually referred to as an overlap of soft wave func- 
tions. It completely neglects effects due to hard PQCD gluon exchanges gen- 
erating the power-law 0((a s /7r) 2 /i 2 ) tail of the nonforward densities at large 
t. Note also that taking the nonforward densities J- a (x, t) with an exponential 
dependence on t, one gets a power-law asymptotics Fj 3oft (t) ~ (— 4A 2 /i)™ +1 
for the Fi(t) form factor, with the power dictated by the (1 — x) n behavior 
of the parton densities for x close to 1. This connection arises because the 
integral ( |162| ) over x is dominated at large t by the region x ~ 4A 2 /|t|. In 
other words, the large-t behavior of Fi(t) in our model is governed by the 
Feynman mechanismEl One should realize, however, that the relevant scale 
4A 2 = 2.8 GeV 2 is rather large. For this reason, when \t\ < 10 GeV 2 , it is 
premature to rely on asymptotic estimates for the soft contribution. Indeed, 
with n = 3.5, the asymptotic estimate is Fj 3oft (i) ~ i~ 4 - 5 , in an apparent con- 
tradiction with the ability of our curve to follow the dipole behavior. The 
resolution of this paradox is very simple: the maxima of nonforward densities 
T a (x,t) for \t\ < 10 GeV 2 are at rather low x-values x < 0.5. Hence, the x- 
integrals producing (t) are not dominated by the x ~ 1 region yet and the 
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asymptotic estimates are not applicable: the functional dependence of Ff oft (t) 
in our model is much more complicated than a simple power of 1/t. 

The fact that the soft overlap model closely reproduces the experimentally 
observed dipole-like behavior of the proton form factor is a clear demonstration, 
that such a behavior does not necessarily reflect the quark counting rules EZtEa 
Ff(t) ~ 1/t 2 valid for the asymptotic behavior of the hard gluon exchange con- 
tributions. Our explanation of the observed magnitude and the ^-dependence 
of Fi(t) by a purely soft contribution is in strong contrast with that of the 
hard PQCD approach to this problem. 



8.3 Wide-angle Compton scattering 

With both photons real, it is not sufficient to have large photon energy to 
ensure short-distance dominance: large-s, small-t region is strongly affected by 
Regge contributions. Hence, having large \t\ > 1 GeV 2 is a necessary condition 
for revealing short-distance dynamics. 

The simplest contributions for the WACS amplitude are given by the s- 
and u-channel handbag diagrams Fig. §b,c. The nonperturbative part in this 
case is given by the proton DDs which determine the i-dependence of the total 
contribution. Just like in the form factor case, the contribution dominating in 
the formal asymptotic limit s, \ t\, u\ — ► oo, is given by diagrams corresponding 
to the pure SD regime, see Fig. ll]a. The hard subgraph then involves two 
hard gluon exchanges which results in a suppression factor (a s /7r) 2 ~ 1/100 
absent in the handbag term. The total contribution- of all two-gluon exchange 
contributions was calculated by_Earrar and Zhangl^l recalculated by Kronfeld 
and Niziclp by VanderhaeghenEiE 2 ] and by Brooks and DixonE3 A sufficiently 
large contribution is only obtained if one uses humpy DAs and 1/fc 2 prop- 
agators with no finite-size effects included. Even with such propagators, the 
WACS amplitude calculated with the asymptotic DA is negligibly small23 com- 
pared to existing data. As argued in Ref. 25, the enhancements generated by 
the humpy DAs should not be taken at face value both for form factors and 
wide-angle Compton scattering amplitudes. For these reasons, we ignore the 
hard contributions to the WACS amplitude as negligibly small. 

Another type of configuration containing hard gluon exchange is shown 
in Fig. Ob. There are also the diagrams with photons coupled to different 



quarks ("cat's ears", Fig. [lip) . Such contributions have both higher order and 
higher twist. This brings in the a s /n factor and an extra 1/s suppression. The 
latter is partially compensated by a slower fall-off of the four-quark DDs with 
t since only one valence quark should change its momentum. 

For simplicity, let us neglect all the suppressed terms and deal only with the 
handbag contributions Fig. ||b,c in which the highly virtual quark propagator 
connecting the photon vertices is convoluted with proton DDs parametrizing 
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Figure 11: Configurations involving double and single gluon exchange. 

the overlap of soft wave functions. Since the basic scale 4A 2 characterizing 
the i-dependence of DDs in our model is 2.8 GeV 2 , while existing data are all 
at momentum transfers t below 5 GeV 2 , we deal with the region where the 
asymptotic estimate (Feynman mechanism) for the overlap contribution is not 
working yet. 

The hard quark propagators for the s and u channel handbag diagrams in 
this case are 



xP + af/2 + Q 
(xP + ar/2 + Q) 2 ~ x~s - (x 2 



xP + af/2 + Q 



x 2 )t/A + x 2 m 2 



(164) 



and 



xP + af/2 - Q 
(xP + ar/2 — Q) 2 xu — {x 2 — a 2 )t/A + x 2 m 2 



xP + af/2 - Q 



(165) 



respectively. We denote s = 2(pq) = s — m 2 and u — —2{pq') = u — m 2 . 
Since DDs are even functions of a, the af terms in the numerators can be 
dropped. Note that it is legitimate to keep 0(m 2 ) and 0(t) terms in the 
denominators: the dependence of hard propagators on target parameters m 2 
and t can be calculated exactly because of the effect analogous to the ^-scaling 
in DISE-3 (see also Ref. 65). Note that the t-correction to hard propagators 
disappears in the large-i limit dominated by the x ~ 1 integration. The t- 
corrections are the largest for a = 0. At this value and for x = 1/2 and t = u, 
corresponding to 90° angle in the center of mass system (cms), the t-term in 
the denominator of the most important second propagator is only 1 /8 of the u 
term. This ratio increases to 1/3 for x = 1/3. However, at nonzero a- values, 
the i-corrections are smaller. Hence, the i-corrections in the denominators of 
hard propagators can produce 10% — 20% effects and should be included in a 
complete analysis. Here, we consider an approximation in which these terms 
are neglected and hard propagators are given by a-independent expressions 
(xP + Q)/xs and (xP + Q)/xu. As a result, the a-integration acts only on the 
DDs f(x, a; t) and converts them into nonforward densities T{x, t). The latter 
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appear through two types of integrals 

f T a {x,t)dx = F?{t) and f T a {x, t) — = B%(t), (166) 

JO JQ x 

and similarly for K, Q,V. The functions F*(t) are the flavor components of the 
usual Fi(t) form factor while R*(t) are the flavor components of a new form 
factor specific to the wide-angle Compton scattering. In the formal asymptotic 
limit \t\ — > oo, the x integrals for Ff(i) and Rx(t) are both dominated by 
the x ~ 1 region: the large-i behavior of these functions is governed by the 
Feynman mechanism and their ratio tends to 1 as \t\ increases (see Fig. ||a). 
However, due to large value of the effective scale 4A 2 = 2.8 GeV 2 , the accessible 
momentum transfers t < 5 GeV 2 are very far from being asymptotic. 

In Fig. [l2|p, the plotted functions are T u (x\t) and T u (x;t)/x at t = —2.5 
GeV 2 . It is clear that the relevant integrals are dominated by rather small x 
values x < 0.4 which results in a strong enhancement of i?"(<) compared to 
F"(f) for \t\ < 5 GeV 2 . Note also that the (p'\ . . . xP . . . \p) matrix elements 
can produce only t as a large variable while (p'\ . . .Q . . .\p) gives s. As a 
result, the enhanced form factors Ri(t) are accompanied by extra s/t factors 
compared to the F®(t) terms. In the cross section, these enhancements are 
squared, i.e., the contributions due to the non-enhanced form factors Ff (t) 
are always accompanied by t 2 /s 2 factors which are smaller than 1/4 for cms 
angles below 90°. Because of double suppression, the terms with F-f (t) can be 
neglected in a simplified approach. 




Figure 12: a) Ratio R% (t) / (t); b) Functions F u {x; t) (solid line) and T u (x; t)/x (dashed 

line) at t = -2.5 GeV 2 . 

The contribution due to the K, functions appears through the flavor com- 
ponents F2 (t) of the F 2 (i) form factor and their enhanced analogues R%(t). 
The major part of contributions due to the /C-type NDs appears in the combi- 
nation R\(t) - (t/4ml)Rl(t). Experimentally, F 2 {t)/F 1 (t) » 1 GeV 2 /|t|. Since 
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R2/F2 ~ Ri/F\ ~ l/(x), R2(t) is similarly suppressed compared to Ri(t), and 
we neglect contributions due to the i?f (t) form factors. We also neglect here 
the terms with another spin-flip distribution V related to the pseudoscalar 
form factor Gp(t) which is dominated by the t-channel pion exchange. Our 
calculations show that the contribution due to the parton helicity sensitive 
densities Q a is suppressed by the factor t 2 /2s 2 compared to that due to the 
T a densities. This factor only reaches 1/8 for the cm angle of 90°, and hence 
the Q a contributions are not very significant numerically. For simplicity, we 
approximate G a (x, t) by T a (x, t). After these approximations, the WACS cross 
section is given by the product 



da 2ira 2 
~dt ~ s 2 



(pq) . GV) 



{pq') (pq) 



R*(t) , (167) 



of the Klein-Nishina cross section (in which we dropped 0(m 2 ) and 0(m . ) 
terms) and the square of the Ri(t) form factor 

Rl(t)=J2 e a[ R l( t )+ R l( t )} ■ ( 168 ) 

a 

In the model of Ref. 25, the effective form factor Ri(t) is given by 
Ri(t) -- 



1 ' e 2 u f v u al (x)+e 2 fr l (x) + 2(e 2 u +e 2 +e 2 ) 







e^l^—. (169) 
x 



The sea distributions are included here assuming that they are all equal 
f sea (x) = fu C ds( x ) = fu d s( x ) an d given by a simplified parametrization 

f sea (x) = 0.5 a;- - 75 (l - x) 7 (170) 

which accurately reproduces the GRV formula for Q 2 ~ 1 GeV 2 . Due to 
suppression of the small- a; region by the exponential exp[xi/4a;A 2 ], the sea 
quark contribution is rather small (~ 10%) even for — t ~ 1 GeV 2 and is 
invisible for —t > 3 GeV 2 . ^ 

Comparison with existing data£3 is shown in Fig. [13]. The curves follow the 
data pattern but are systematically lower by a factor of 2, with disagreement 
becoming more pronounced as the scattering angle increases. Since several 
terms were neglected each capable of producing up to a 20% correction in 
the amplitude, the agreement between curves and the data may be treated as 
encouraging. The most important corrections which should be included in a 
more detailed investigation are the t-corrections in the denominators of hard 
propagators and contributions due to the "non- leading" K,,Q,V nonforward 
densities. The latter, as noted above, are usually accompanied by t/s and t/u 
factors, i.e., their contribution becomes more significant at larger angles. The 
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Figure 13: Angular dependence of the combination s 6 (da/di). 



^-correction in the most important hard propagator term l/[xu— (x 2 — a 2 )t/A+ 
x 2 m,p} also enhances the amplitude at large angles. 

Note that the curves for the combination s 6 (da/dt) taken for the initial 
photon energies 2, 3, 4, 5 and 6 GeV intersect each other at cm ~ 60°. This 
is in good agreement with experimental data of Ref. 53 where the differential 
cross section at fixed cms angles was fitted by powers of s: da/dt ~ s - "^) with 
n cx P(60°) = 5.9±0.3. The curves of Ref. 25 correspond to n soft (60°) w 6.1 and 
n soft (90°) w 6.7 which also agrees with the experimental result n exp (90°) = 
7.1 ±0.4. 

This can be compared with the scaling behavior of the asymptotic hard 
contribution: modulo logarithms contained in the a s factors, they have a uni- 
versal angle-independent power n hard (#) = 6. For cms = 105°, the experimen- 
tal result based on just two data points is n exp (105°) = 6.2 ± 1.4, while our 
model gives n soft (105°) w 7.0. Clearly, better data are needed to draw any 
conclusions here. 
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9 Concluding remarks 

In this paper, I described the basic elements of the theory of generalized parton 
distributions. For uniformity of presentation, I heavily relied on the approach 
developed in my papers, see Refs. 16-18 and 21-25. In this concluding section, 
I briefly list other developments in the theory of GPDs and its applications, 
not covered in the present review. Additional references can be found in the 
reviews and review sections of the original papers cited in Refs. 11, 14, 17, 
19-21, 45, 66, 50, and 67. 

Introduction of GPDs and factorization. In various ways (and under dif- 
ferent names: oflkEorward., nonforward, nondiagonal, off-diagonal, etc.) GPDs 



iiwawxLjac 

were introduced tSliEllla'ESEJ as nonperturbative functions describing the non- 
perturbative part of the factorized representation for the amplitudes of hard 
elastic electroproduction processes. The PQCD factorization for DVCS and 
meson production processes was discussed in Refs. 17, 20, 68, 69. One loop 
corrections to the DVCS amplitude were calculated in Refs. 70-72. 

Evolution of G-FPfi .Evolution equations for GPDs obtained originally at 



one loop leveltlrllZIOEJ'EiJ were used fip£.numerical studies of ev^piflaip using the 

;sE-30 ; 



orthogonal polynomials techniquesEJO and direct integrationlI!3Qi2j Two-loop 
evolution was investigated in Refs. 78-80. 

Applications of GPDs to hard electroproduction processes. There is a grow- 
ing literature devoted to practical aspects of using GPDs in the description 
of hard electroproduction processes. In particular, deeply virtual Compton 
scattering is discussed in Refs. 11, 44, 66, 62, 81-84. The PQCD approach to 
meson electroproduction at large Q 2 formulated. p. Refs. 20, 16, 17 was anplied. 
to particular channels: p-meson productionpJ'L^rEB, pion production(3E3E3 
and, finally, electroproduction processes accompanied by excitation of decu- 
plet baryonsa 
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